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ABSTRACT 


Young open clusters located in the outer Galaxy provide us with an oppor¬ 
tunity to study star formation activity in a different environment from the solar 
neighborhood. We present a UBVI and Ha photometric study of the young 
open clusters NGC 1624 and NGG 1931 that are situated toward the Galactic 
anticenter. Various photometric diagrams are used to select the members of the 
clusters and to determine the fundamental parameters. NGG 1624 and NGG 
1931 are, on average, reddened by {E{B — V)) = 0.92 ± 0.05 and 0.74 ± 0.17 
mag, respectively. The properties of the reddening toward NGG 1931 indicate 
an abnormal reddening law (Ryd = 5.2 ± 0.3). Using the zero-age main se¬ 
quence fitting method we conhrm that NGG 1624 is 6.0 ± 0.6 kpc away from 
the Sun, whereas NGG 1931 is at a distance of 2.3 ± 0.2 kpc. The results from 
isochrone htting in the Hertzsprung-Russell diagram indicate the ages of NGG 
1624 and NGG 1931 to be less than 4 Myr and 1.5 - 2.0 Myr, respectively. We 
derived the initial mass function (IMF) of the clusters. The slope of the IMF 
(rNGC 1624 = —2.0 ± 0.2 and Fngc 1931 = —2.0 ± 0.1) appears to be steeper than 
that of the Salpeter/Kroupa IMF. We discuss the implication of the derived IMF 
based on simple Monte-Garlo simulations and conclude that the property of star 
formation in the clusters seems not to be far different from that in the solar 
neighborhood. 

Subject headings: open clusters and associations: individual (NGC 1624, NGC 1931) 
~ local interstellar matter - stars: luminosity function, mass function 
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INTRODUCTION 


Young open clusters are useful objects to study the star formation process because 
about 80 percent of th e stars in star forming regions (SFRs) are found in clusters with 


at least 100 members ( Lada fc Lad 


2003 


Porras et al 


20031 ). The relation between star 


formation activity and environmental condi 

;ions is one 

in star formation studies ( 

Garamazza et al. 

00 

o 

o 

CM 

CM 

1 

o 

CM 


20121) . A basic diagnostic tool for 


understanding star formation p rocesses 


concept was hrst introduced by 


is the stellar initial mass function (IMF). The 


Salpeterl fjl955l) . If the stellar IMF has a universal shape 


in all SFRs, there must be a fundamental factor controlling the star formation process. 
A number of effort s have been devoted to conhrming the universality (see review of 


Bastian et al. 


20101) . however, there is as yet no firm conclusion concerning the universality 


or diversity of the shape. In this context, the IMF of young open clusters formed in various 
star-forming environments may give a clue to the dependence of star formation processes 
on environmental conditions. As part of an attempt to study this issue, we investigated the 
young open clusters NGC 1624 and NGC 1931 in the outer Galaxy, which is probably a 
low- metallicity environ ment, according to the gradient of metal abundance in the Galactic 


disk (lYong et al. 


29121 ). 


The young open cluster NGG 1624 is surrounded by the H ii region Sh2-212 with 
a shell like structure. The cluster comprises several high-mass mai n sequence (MS 


stars, as well as a large number of pre-main sequ ence (PMS) stars (IJose et al. 


The most luminous star is NGG 1624-2 (07f?p 


Walborn et al. 


29191: 


2911). 


Sot a et al 


291ll) 


located near the cluster center. This star with strong X-ray emission is known to be a 


slowly rotating star with a very strong magnetic field flWade et al. 


probable PMS mer nbers one-fifth ha v e circ umstellar disks (iJose et al. 


radio observations. 


Deharveng et al. 


de et al. 

29 

12 

). c 

(Jose et al. 

2911) 


(129981) found an ultra-compact H ii (UGHII) region 
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at the border of Sh2-212. A high-mass PMS star (No. 228 of Deharveng et ah 2008 ) 


seems to be associated with the UCHII. A 
the fundamental parameters 


20061 : 


Deharveng et ah 


2008 


Moffat et al 


Jose et al 


ew photometric studies o: 


1979 


Chini fc Wink 


1984 


NGC 1624 provid ed 


Suiatha fc Babu 


201 ll) . The reddening of the cluster is about 


E{B — V) = 0.7 - 0.9 r nag. The previous l y determined distance 


the range 6.0 - 6.5 kpc 


Jose et al 


2011 ). while 


(Moffat et al. 


1979; 


IMF of the cluster was investigated by 


Suiatha fc Babu 


to the cluster was in 


2006 


Deharveng et al. 


2008 


Chini fc WinkI ( 19841) obtained a far dif feren t value of 


Suiatha &: Babul f 20061) and 


its s 


2001 


ope wa s reasonably consistent with the Salpeter/Kroupa IMF flSalpeter 


Jose et a 


0 .3 kp c. The 
J2OI1I). and 


1955 


Kroupa 


2 OO 2 I ). It implies that although NGC 1624 is located in the outer Galaxy, the nature 


of the star formation activity is similar to that found in the solar neighbourhood. 

The nearer and younger open cluster NGC 1931 is associated with the glowing 


nebula Sh2-237. A 


fjPandev et al. 


ew bright MS stars and many PMS stars constitute the cluster 


2013al). The main io nizing source of t 


the BO.5 star (IGlushkov et al. 


le bright nebula is 


PMS star candidates were identified b y 


Sky Survey (2MASS; 


Skrutskie et al. 


bought to be 


';wo B2 stars 

Pandev et al. 

2013a 

). A number of 

Pandev et al. 

(2013a 

) from the Two Micron All 


200611 catalog and Spitzer InfraRed Array Camera 


(IRAC) photometry. The spatial distribution of members seems to be divided into tw o 


2009 


Pandev et al. 


stellar groups, a northern and a southern group fiBonatto fc Bica 
2013al) . The fundarnental parameters of the cluster were obtained from photometric 


2009 


studies flMoffat et a 


Pandev et al 


1979; 


Pandey fc Mahra 


1986 : 


Bhatt et al. 


1994 : 


Bonatto fc Bica 


2013al) . According to these studies, the reddening and distance of 


the clu ster are EjB 


study fjPandev et al. 


V) = 0.5 - 1.0 mag and d = 1.8 - 2.4 kpc, respectively. A recent 


2013al) investigated the properties of dust toward the cluster using 
polarimetric and p hotometric data. The reddening law toward the cluster estimated from 


the Serkowski law fjSerkowski et al 


I 975 I) and color excess ratios indicated that the size 
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distribution of dust grains may be different from that of the general diffuse interstellar 
medium (ISM). They also derived the IMF of the northern, southern, and the entire 
cluster regions, respectively. The slope of the IMF appears to be shallower than the 
Salpeter/Kroupa IMF. 


The present work on NGC 1624 and NGC 1931 is the hfth paper of the Sejong Open 
cluster Survey (SOS) project which was initiated to provide homogeneo us photomet r ic data 


Sung et al.l (2013a) 


for many open clusters. The overview of the SOS project can be found in 
(hereafter Paper 0). Gomprehensive studies of several open clusters NGG 235 3, IG 1848, 


and NGG 1893 were published as part of the project flhim et ah 


2011 


2014alJb h In this 


work, we revisit the reddening law as well as the fundamental parameters of these clusters 
in a homogeneous manner. The IMF of the clusters is also studied in order to investigate 
the property of star formation activity in the outer Galaxy. The observation and reliability 
of our photometry are described in Section 2. In Section 3, we present several fundamental 
parameters of the clusters obtained from photometric diagrams. The reddening law toward 
two clusters is also discussed in this section. We construct the Hertzsprung-Russell diagram 
(HRD) in Section 4 and derive the IMF of the clusters in Section 5. Several discussions on 
the age spread of PMS stars are made in Section 6. Finally, the comprehensive results from 
this study are summarized in Section 7. 


2. OBSERVATION 


2.1. AZT-22 1.5m Telescope of Maidanak Astronomical Observatory 


The observations of NGG 1624 was made on 2006 November 24, using the AZT-22 
1.5m telescope (f/7.74) at Maidanak Astronomical Observatory (MAO) ii i Uzbekistan. A ll 


imaging data were acquired using Fairchild 486 GGD camera (SNUGam; 


Im et ah 


2010h 
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with the standard Bessell UBVI filters (Bessel] 1990 ). The field of view 


18(1 X 18(1. The characteristics of the CCD chip have been described by 


fFOV) is a 


Dont 


him et al. 


fl2nn8h 


in detail. The observations comprised 8 frames that were taken in two sets of exposnre 
times for each band - 3 and 60s in J, 3 and 180s in V, 5 and 300s in B, and 15 and 600s 
in U. To transform instrumental magnitudes to s tandard magnitudes and colors we made 


observations of several equatorial standard stars (IMenzies et al. 


199 if) over a wide air mass 


range of X ~ 1.2 - 2.0 on the same night. The observations are summarized in Table [H 
and the left-hand panel in Figure [1] displays the hnder chart (dashed line) for the observed 
stars brighter than D = 18 mag. 


Al 

n 


pre-processing to remove instrume ntal artifac t s was carried out using the 


IRAFtl/CCDRED packages as described in 


him et al. 


fj2008l) . Simple aperture photometry 


was performed for the standard stars with an aperture size of 14.0 arcsec (26.3 pixels). 
Atmospheric extinction coefficients were determined from the photometric data of the 
standard stars with a weighted least-square method. The coefficients are presented with 
photometric zero points in Table [2l Using IRAF/DAOPHOT we conducted point spread 
function (PSF) photometry on stars in the target images using a small fitting radius of 1 
full width at half-maximum (FWHM; <1.0 arcsec), and then aperture correction was done 
using the aperture photometry of bright, isolated stars with a photometric error smaller 
than 0.01 mag in individual target images. The instrumental magnitude was transformed 
to the standard magnitude and colors using the transformation equations below (Paper 0): 


Mx = mx- {hx - k2xC) ■ X + rjx ■ C + ax ■ UT + (1) 

^ Image Reduction and Analysis Facility is developed and distributed by the National 
Optical Astronomy Observatories, which is operated by the Association of Universities for 
Research in Astronomy under operative agreement with the National Science Foundation. 
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where Mx, m\, ki\, k 2 \, rjx, C, X, a\, UT, and are the standard magnitnde, instrumental 
magnitude, the primary extinction coefficient, the secondary extinction coefficient, 
transformation coefficient, relevant color, air mass, time-variation coefficient, time difference 
relative to midnight, and photometric zero point, respectively. We adopted the SNUCa m 


trans; 


(120091). 


brmation coefficients {tjx), which were recently refined from those in 


Tam et ah 


2.2. Kuiper 61” Telescope of Steward Observatory 

The observations of NGC 1624 and NGC 1931 were made on 2011 October 29, using 
the Kuiper 61” telescope (f/13.5) of Steward Observatory (SO) on Mt. Bigelow in Arizona, 
USA. We obtained images with the Montlk GGD camera and 5 hlters (Bessell U, Harris 
HU, Arizona /, and Ha). The FOV is about 9'.7 x 917. Later, additional images of NGG 
1931 were taken on 2013 February 5, without the Ha hlter. Since typical seeing was I'.'O - 
2'.'0 (7 ~ 14 pixels), the 3x3 binning mode (0742 per pixel) allowed us not only to obtain 
an appropriate FWHM (~ 3 pixels) for PSF photometry but also to decrease the readout 
time (~ 10s). In order to obtain reliable transformation relations for this photometric 


Kilkennv et ah 

1998 

-V—QV-., ..X.. V.. X.. 

) as well as a large number of equatorial standard stars ( 

Menzies et ah 


199ll) at air masses of 1.2 - 2 during the run. In ad dition, a few Landolt standard star 


fields (Rubin 149, 152, and PG 2213 


Landolt 


19921), which contain both extreme red and 


blue stars within a small region, were also observed to obtain the secondary extinction 
coefficients. The observations are summarized in Table [H and the finder chart (solid line) 
for the stars brighter than U = 18 mag is shown in Figured] 


The Mont4k photometric system of the Kuiper 61” telescope was used in this survey 
project for the hrst time. Understanding the characteristics of the photometric system 

















is essential for obtaining reliable photometric data. We address the characteristics and 
transformation relations for the Mont4k CCD photometric system in the Appendix. 
Atmospheric extinction coefficients and photometric zero points obtained from the 
apertnre photometry of the standard stars are presented in Table [2l PSF photometry and 
aperture correction for the images of NGC 1624 and NGC 1931 were carried out using 
the same procedures as described in the previous section. The instrumental magnitude 
was transformed to the standard magnitude and colors using the transformation equations 
addressed in the Appendix. 


2.3. The Consistency and Completeness of the Photometric Data 


We have conhrmed from a series of studies flLim et al 


2011 


2014a b) that the 


photometric data obtained from the photometric system of the AZT-22 1.5m telescope 
at MAO are well-tied to the Johnson-Cousins standard system. But we need to check 
the homogeneity of the photometric data obtained with the Mont4k CCD camera of the 
Kuiper 61” telescope. Since NGC 1624 was observed with the two different photometric 
systems, it was possible for us to compare our photometric data directly. Figure H] shows 
the comparison between a few sets of modern CCD photometry. Bold dots (black) in the 
upper panels represent the differences between the photometric data obtained at MAO and 
SO. We confirmed that our photometry is in good agreement within 0.01 mag (see also 
Table E]). The imaging data for NGC 1931 were taken with the Mont4k CCD camera of the 
Kuiper 61” telescope at two different epochs. The internal consistency of the two data sets 
was examined as shown by bold dots in the lower panels of Figure [2l The photometry is 
well consistent with each other wi thin 0.01 mag. We averaged the data sets for each cluster 


using a weighted average scheme fISung &: Lee 


199, 5h 


A couple of photometric studies with a modern CCD camera have been made for 
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NGC 1624 (e.g. ISuiatha &: Babul I 2 OO 6 I : I Jose et alJboill) . The combined photometric data 
for the cluster were compared with those of the previous studies as marked by squares 


and triangles in t 


le up per panels of Figure [2] (see also Table E]). The photometric data of 


Suiatha fc Babul (120061 1 are far different from ours, particularly the V magnitudes. All their 


colors are systematically redder than the others. There is also a color-dependent trend for 
U-B< 0.7. 


The difference between the photometry of 


Jose et ah 


(1201 ll) and ours is acceptable 


in V and B — V, whereas the other colors show systematic differences. Although their 
U — B color is in a good agreement with ours for U — B < 0.4, for redder stars their U — B 
colors are systematically redder than ours. Their V — I color appears bluer than that of 
the others. We note that the previous sets of CCD data for NCC 1624 are likely to involve 
problem s in the photometric s tandardization. On the other hand, the photometry of NCC 


1931 by 


Pandev et ah 


fl2013all is in good agreement with ours. We present a summary of 


the comparison in Table El 


Statistical analysis based on photometric data requires the completeness of the 
photometry to be at least 90 percent. In the case of sparse open clusters, the completeness 
limit is not very different from that obtained in nearby field region, and that relates to 
the turnover magnitude in the luminosity function of the observed stars. We estimated 
the completeness of our photometry from the turnover in the luminosity function of all 
the observed stars, by assuming a linear slope across the entire magnitude range. Our 
photometry is complete down to C = 19 mag (~ 2 Mq) for NCC 1624 and 20 mag 
(~ 1 Mq) for NCC 1931. However, the completeness limit may be an upper limit because 
the bright nebulae surrounding the young open clusters may prevent the detection of faint 
stars. The photometric data from this work are available in the electronic tables (Table 0] 
andE]) or from the authors (BL or HS). 
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3. PHOTOMETRIC DIAGRAMS 

Photometric studies for young open clusters are based on the two-color diagram 
(TCD) and color-magnitude diagrams (CMDs). With well-calibrated empirical relations 
and stellar evolution models, fundamental parameters, such as reddening, distance, and 
age, can be determined from these photometric diagrams. The fundamental parameters 
allow us to investigate the local spiral arm structure in the Galaxy (Paper 0) as well as 


the properties and evolution of the dust grain 


2003: 

him et al. 

2011 , 

2014a 

b; 


Sung fc Bessell 


hxough the reddening law fjPandev et ah 


2014 Paper 0, etc). In addition, the age 


distrib ution of stars can be used to study the star form a tion history within a s 


region flPalla fc Stabler 


1999 


2ni4b 


Sicilia-Aguilar et al. 

2004; 

Sung & Bessell 

2010 

him et al. 


, etc) and the pattern speed of the Galactic spiral arms flDias fc Lepine 


;ar-forming 


200511 . These 


studies should be based on reliable membership selection. In the following sections, we 
present membership selection criteria, the reddening law, and fundamental parameters of 
NGG 1624 and NGG 1931 with TGDs and GMDs. 

3.1. The Extent of NGC 1624 and NGG 1931 


As most open clusters are unbound stellar systems, i.e. expanding systems 


(iLada fc Lada 


20031 1. it is believed from theoretical approaches, that gas expu 


affects the dynamical evo lution of the remaining clusters signihcantly flTutukov 


Goodwin &: Bastian 


Sion 


1978 : 


20061) . In some SFRs, new generation stars are still forming in the 


natal clouds swept away 


(e.g. 


Sicilia-Aguilar et al 


3 V the stellar wind and the UV radiation fro m high-mass stars 


2004; 

Koenig et al. 

2008; 

him et al. 

to 

o 

g _ 


structure is basically a projection onto the sky. Hence, it is difficult to determine the 
physically meaningful extent of young open clusters. Alternatively, radial surface density 
prohle (RSDP) gives an area encompassing almost all the members of a cluster. This is a 
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fundamental method to define the size of clusters. 


We determined the radius of NGC 1624 assuming a circular shape. The center of the 
cluster was set to the position of the brightest star NGC 1624-2. The stars observed in 
the V band were counted within concentric rings increased by 0^5. The RSDP of NGC 
1624 is plotted in Figure O The error in the surface density was assumed to follow Poisson 
statistics. The surface density decreases with the increasing distance from the center until 
it indistinguishably converges to the surface density of held stars. The distance from the 
center to the convergent point (~ 2!5) was assigned as the radius of the cluster. 


Jose et ah 


(1201 ll) also obtained a consistent value (2f0) with ours using their optical and 2MASS 
near-infrared (NIR) data. 


The extent of NGC 1931 could not be determined from the RSDP. 


Pandev et ah 


fj2013al) presented the spatial distribution of PMS members in their hgure 25. The 


spatial distribution could not be conhned to a small specihc region within our small FOV 
(~ 10' X 10') because the PMS stars scattered across the emission nebula fSh2-237). 
Given that no PMS member has been identihed northwards from 6 ~ 34° 18' fjPandev et ah 


2013aJ), we assigned the southern region covered by our observation as likely to encompass 


all the members of the cluster. The area of these clusters is used to normalize the IMF as 
well as to isolate members. 


3.2. Membership Selection 

The early-type stars (O - B-type) can be unambiguously identified in the (1/ — R, 

B — V) TCD as they have very blue U — B colors (see hgure 11 of Paper 0). In addition, 
the intrinsic colors and absolute magnitude of these stars have been well calibrated in 
the optical passbands, so that the reddening and distance of the stars can be reliably 
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determined by comparing the observed colors and magnitude with the intrinsic relations. 
These procedures are described in detail in the following sections. We identihed early-type 
members within specihc reddening and distance ranges based on the {U — B, B — V) TCD 
(Figure Hj) and CMDs (Figure [5]). The membership selection criteria for the early-type 
main sequence (MS) stars of NGC 1624 are (1) Id < 18 mag, 0.5 < B — V < 0.9, 

-1.0 <U -B < 0.5, and Q' < -0.4, where Q’ = {U - B) - 0.72(S -V)- 0.02bE{B - Vf 
(Paper 0), (2) E{B — V) > 0.8 mag, (3) for late-B type members (—0.4 < Q' < —0.2) 
the individual distance modulus should be between (Vq — Mv)c\ — 0.75 — 2avo-Mv and 
(Vo — My)c\ + 2avQ-Mv-i where (Vq — Mv)c\ and av^-Mv are the distance modulus of NGC 
1624 and the standard deviation of the distance modulus, respectively. The factor -0.75 is 
introduced to take into account the effect of equal mass binaries. The FOV of SNUCam 
(1811 X 1811) is wide enough to include the distributed population, most of which may not 
be associated with NGC 1624. In order to prevent the inclusion of such stars we isolated 
only stars within a radius of 215 from the brightest star NGC 1624-2. The star ID 1609 
{V = 14.94, V — I = 1.05, B — V = 0.87, U — B = 0.24 ) was rejected from the member 
list because the star is too bright to be the member of the cluster at a given color in the 
CMDs. These processes of the membership selection for early-type MS stars were carried 
out iteratively. A total of 14 stars were selected as the early-type MS members of NGC 
1624. 

Similar criteria were used to select the early-type MS members of NGC 1931. The 
criteria are (1) V <17 mag, 0.2 < B — V < 0.9, —1.0 < U — B < 0.4, and Q' < —0.4, (2) 
E{B — V)> 0.5 mag, (3) 10.8 mag < {Vq — My)* < 12.3 mag for late-B type members. 
Several held interlopers were identihed from these criteria. We checked the color excess 
ratio of the most probable members {U — B < —0.1) and excluded outliers with signihcantly 
diherent color excess ratios from those of the most probable members. The abnormal color 
excess ratios of these outliers are caused by the fact that less reddened foreground F or 
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G-type stars would be regarded as highly reddened early-type stars from the color-cuts. 
The star ID 314 (V = 13.11, V — I = 0.53, B — V = 0.36, U — B = —0.31) was originally 
selected as an early-type member, however, its spatial position {Aa ~ —3'.23, A6 ~ -I-5D1) 
is outside the cluster boundary. We regarded the star as a mid-B-type field interloper at 
a similar distance, and therefore the star was excluded from the member list. A total of 
14 stars were assigned as early-type MS members of NGC 1931. It is worth noting that 
BD-l-34 1074, which is one of the two brightest stars in NGG 1931, was resolved into four 
stars, of which three (ID 857, 865, and 872) are also the early-type members identihed here. 


Most low-mass PMS stars have a circumste 


significant accretion a ctivity flLim et ah 


model for PMS stars (jUchida fc Shibata 


2014a 


1985 


lar disk, and some of these disks exhibit 


). According to t 


Bertout et ah 


le standard accretion 


1988 


Konig 


19911), material 


channeled from the circumstellar disk falls onto the surface of the central star along its 
magnetosphere. Various emission lines and hot continuum excess emi ssion arise from 


accre t ion columns, pre shock infall region, and the heated photosphere flGalvet fc Gullbring 


19981 : 


Hartmann 


19991 ). The generated energy is mainly released at ultra violet (UV) 


wave l engths, so that the accreti on luminosity is related to the U magnitude flGullbring et ah 


19981: 


Galvet fc Gullbring 


199811 . Hence, a few PMS stars with strong UV excess emission 


can be identified in the U band photometry. There is a dotted line parallel to the reddening 
vector in the {U — B, B — V) TGD (the left-hand panels of Figure H]). The stars with a bluer 
U — B color than the dotted line at a given B — V color are either early-type stars or UV 
excess stars. The criteria for the PMS stars with UV excess emission are (1) e{U — B) < 0.1, 


(2) B — U > 1 anc 


locus (jSung et ah 


U — B < 0.72[(i? — U) -1- 0.5)] — 0.9, (3) stars within the empirical PMS 


200811 in the (U, V — I) GMD (see Figure [5]1. The last criterion prevents 


the inclusion of background early-type stars. We found 3 and 4 PMS stars with UV excess 
emission in NGG 1624 and NGG 1931, respectively. 
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Ha photo metry provi d es a 


3 Myr). Since 


Sung et al. 


good criterion to find PMS stars in young open clusters (< 


fjl997l) achieyed the successful detection of many PMS stars in 


of studies (e.g. 

Sung et al. 

1998, 

2000 ; 

Park et al. 

2000 


Park & 

Sung 

2002 

i: Sung 

et al. 

20Q4; 

Sung & Bessell 

20Q4; 

Sung et al. 

2008: 

Sung et al. 

2013h: 

him et al. 

2014a 

d). In 


order to id entify Ha emissio n stars photometrically we haye defined the Ha index as Ha 


{V + /)/2 fISung et ah 


2000|). As shown in the right-hand panels of Figure IU stars with -0.2 


(dashed line) or -0.1 mag (dotted line) smaller Ha index than the empirical photospheric 
leyel (solid line) of normal MS stars are selected as Ha emission stars or candidates. The 
Ha emission stars and candidates ID 1897 {V = 18.68, B — V = 1.17, U — B = 0.41), 

ID 2326 {V = 18.97, B -V = 1.62, U - B = 1.42), ID 2356 {V = 19.11, B -V = 1.36, 

D - H = 1.14), and ID 2588 (V = 19.12, B - V = 1.17, U - B = 0.44) in the FOV of 
NGC 1624 were regarded as foreground stars because their colors are similar to those of 
less reddened foreground stars. In addition, the radius (215) of the cluster was also used to 
isolate the members. We found 9 Ha emission stars and candidates associated with NGC 
1624. In the same way, we found 28 Ha emission stars and candidates in NGC 1931, one of 
which (ID 596; V = 19.30, B — V = 1.36, U — B = 1.00) is likely a foreground late-type star. 

The dust emission from circumstellar disks of PMS stars is prominent at infrared (IR) 
wavelength, particularly the mid-infrared (MIR). A number of young stellar objects in 


various star-forming regions have been identified t 
the Syitzer space telesc ope fjCaramazza et al. 


i rough extensive imaging surveys with 


2008 


Gntermnth et al. 


2008 


Sung et al. 


2008 : 


Koenig et al. 


20091. etc). We used the Spitzer Galacti c Legacy Infrared Iv 


Survey Extraordinaire 360-degree catalog (GLIMPSE360; 


Whitney et al. 


2008 


id-Plane 


201 ih to 


identify PMS members with MIR excess emission. The GLIMPSE360 suryey is a “Warm 
Mission”, so that only the 3.6 and 4.5 /rm bands were available. Unfortunately, NGC 1624 
{I = 155?356, b = 2?616) was not covered in the survey program because of its high Galactic 
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latitude. The catalog allows us to select PMS members with a circumstellar disk in NGC 
1931. A total of 1249 optical counterparts in the catalog were found within a matching 
radius of 1"0. We attempted to identify the PMS members in the ([3.6] — [4.5], Id — I) 
TCD as shown in Figure [6l The majority of the Ha emission stars exhibited MIR excess 


emission. Stars with a [3.6] — [4.5] color larger than 0.2 mag (da shed line) were se 


as PMS member candidates, and then the empirical PMS locus fjSung et ah 


2008 


ected 


was 


used to isolate the PMS members among the candidates. From this procedure, a total of 
54 PMS stars with MIR excess emission were selected as members of NGC 1931. Although 
it is difficult to make a detailed classihcation due to the absence of information in other 
IRAC bands (5.8 ari d 8.0 /rm), most of them may be Class II objects (see also Figure 18 of 


Pandev et ah 


2013ah . 


It is a well known observational fact that PMS stars, particularly classical T-Tauri 

ects with an amplitude of 0.1 - 2 mag in the form of irregular 


(CTTS) stars, are variab 

e obi( 

variation ( 

Grankin et ah 

2007) 


20071) . The observations of NGC 1624 and NGC 1931 were made 


at two different epochs. The time difference between the hrst and the second observation 
was about ~ 4.9 and 1.2 years for NGC 1624 and NGC 1931, respectively. Since the 
photo metric errors froni the averaging process used in section 2.3 (originally from equation 


(2) of 


Sung fc Lee 


1995h represent the consistency of magnitude and colors from several 


observations, a photometric error larger than that expected from the distribution of errors 
with magnitude may be attributed to a genuine variation in brightness. We investigated 
the photometric error of the stars observed more than twice as shown in Figure [71 The 
photometric errors {(Tv) in a given magnitude bin (AR = 1 mag) were averaged, and then 
its standard deviation {as.n.) was used as a criterion for the variability of individual stars. 
If stars have variations larger than 0.03 mag in brightness and ay > {<tv) + 3<ts.d. in the 
range of R = 10 - 21 mag, we assigned the stars as variables. The star ID 727 in NGC 1931 
showed a very large variation (R = 16.34 mag on 2011 Oct. 29 and R = 14.86 mag on 2013 
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Feb. 5). 

We found 4 variable stars within the radius of NGC 1624. The variable star ID 1733 
has a red U — B color similar to that of foreground stars. The U — B color of the star in 
2011 was bluer by 0.24 mag than that observed in 2006. The star may be an active late-type 
star in front of the cluster. The V magnitude and V — I color of the other variable stars 
appear to be commensurate with those of the Ha emission stars in the PMS locus. A total 
of 16 stars were identified as variable stars in NGC 1931. The majority of the variable stars 
were crowded into the cluster center. Their V magnitudes and colors are similar to those of 
other PMS members with Ha and MIR excess emission, some of which (~ 44 percent) are 
indeed UV excess, Ha emission, or MIR excess emission stars. Thus, these variable stars 
were assigned as PMS members of NGC 1931. 


A total of 28 stars (14 early-type and 14 PMS stars) in NGC 1624 and 85 stars 


(14 early-type anc 
1773 (No. 228 of 


71 PMS stars) in NG C 1931 were selected as members. The star ID 


Deharveng et ah 


20081 ) in NGC 1624 is the only star selected using 2 


membership selection criteria (UV excess and Ha emission), and we found 24 out of 71 
PMS members in NGC 1931 satisfied more than 2 selection criteria. The membership 
selection for early-type MS members is likely complete, however only a small fraction of 
the PMS stars may have been selecte d as members of t he clusters. According to a study 


of the young open cluster NGC 1893 flLim et al 


3 , 


2014m, the detection efficiencies of Ha 


photometry, Spitzer MIR, and Chandra X-ray data turned out to be about 10, 24, and 85 
percent for PMS stars (> 1 Mq), respectively. This implies that X-ray emission from PMS 
stars is the most efficient criterion to identify the remaining PMS stars. Although NGC 
1624 has been observed using the Chandra X-ray Observatory (ObsID 7473, PI Garmire), 
the exposure times were not long enou gh to detect PMS members. Using the published 


X-ray source list of 


Evans et al. 


fl2010[) . the optical counterparts of X-ray emission sources 
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and candidates were searched for with matching radii of 1.0 and 1.5 arcsec, respectively. 
We fonnd only 2 X-ray emission sonrces and 1 candidate with optical connt erparts. The 


brightest early-type me mber NGC 1624-2 is known to be an X-ray emitter flEvans et ah 


2010 


Wade et ah 


201211 . and the others tnrned ont to be held interlopers given their colors. 


On the other hand, X-ray observations for NGC 1931 have not yet been made. Extensive 
X-ray imaging observations are clearly reqnired in order to stndy the natnre of PMS stars 
in detail based on complete membership lists. As of now, we merely anticipate that a few 
hnndreds of members may exist in each cluster. A discussion on the approximate number is 
addressed in a later section. 


3.3. Structure of the Clusters 


The structure of young open clusters gives us clues to dynamical evolution in the 


early stages of c^ 


(lElmegreen et al 


uster formation as well as the formation process of stellar clusters 


200011 . Many efforts have been made to study the structure of open 


clust ers using the 


techniques of surface density distribution and minimum spanning 


tree ( 

Sung 

et al. 

1999; 

Sharma et al. 

2007 

; I 

<oenig et al. 

2008; 

Cute 

rmnth et al. 

2009 

• 

Kook et al. 

2010 ; 

Jose et al. 

2011 ; 

Sung et al. 

2013b: 

him et al. 

2013; 

Pandev et al 

. 2013a 


etc). In this section, we describe the apparent structure of NGC 1624 and NGC 1931 using 
the surface density distribution, with additional information from previous studies. 

The appearance of the HII region Sh2-212 associated with NGC 1624 is close to a 
symmetric shell structure. An inner region hlled with hot i onized gas and a shocked outer 


region constitute the apparent gas structure (see hgure 1 of 


Jose et al 


1201 ll ). A molecular 

hlament with different velocity components in the range of -32 to -37 km s~^ surrounds 


the southern part of the cluster, stretching out to the north-west flDeharveng et al. 


200811 


The hlament incubates at least hve clumps, one of which is likely to be associated with 
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an UGHII in t 

le western part of the cluster. The brightest star NGG 1624-2 (07f?p - 

Walborn et al. 

2010 ; 

Sota et al. 

2011 

) is located near the center of the cluster. Other 


members are distributed in the vicinity of the star as shown in the left-hand panel of 
Figure m We obtained the surface density distribution of NGC 1624 using stars observed in 
the V band (contour). The concentration of stars appears high in the center of the cluster, 
and the spatial distribution approximates to a circular shape. 


of 


On the other hand , NGC 1931 is enclosed within a dusty molecular cloud (see hgure 1 


Pandev et al. 


2013al) . The natal cloud also exhibits a hierarchical structure. The inner 


region is hlled with hot ionized gas produced by the radiation from a few bright stars while 
polycyclic aromatic hydrocarbon molecules are glowing in the outer region. The right-hand 
panel in Figure |8] shows the spatial distribution of cluster members. The early-type MS 
members appear to be divided into two groups, a northern and a southern group. A high 
stellar density region is seen between the groups, implying that mass segregation among the 
members may not yet be established. While many variable PMS members are located in 
the dense region, the majority of Ha emission stars are found in the vicinity of the southern 
group. The PMS members with MIR excess emission are spread out across the whole 
region, however the stars in the northern part of the cluster show a weak concentration. 
The surface density distribution (contour) refle cts such an elongated shape. The apparent 


shape is in a good agreement with the result of 


Pandev et al 


fj2013al) . These observational 


properties may be related to the star formation history within NGG 1931. 


3.4. Reddening and the Reddening Law 

As most open clusters are formed in the Galactic plane, a region where interstellar 
matter is unevenly distributed, reddening corrections are crucial to obtaining reliable 
physical quantities. The interstellar reddening is basically determined by comparing 
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an observed color with the intrinsic one. The {U — B, B — V) TCD is a very 
nsefnl tool because the reddening vector has been well established in the diagram 
[E{U — B)/E{B — V)= 0.72 + 0.02577(5 — V)- Paper 0]. The intrinsic color relations of 
Paper 0 were adopted to obtain the reddening of the individual early-type members (see 
table 1 in Paper 0). 


The reddening of NGC 1624 determined from 14 early-type members is in the range 
of E{B — V) = 0.83 - 1.01 mag, and the mean value is {E{B — V)) = 0.92 ± 0.05 (s.d.) 
mag. The re sult is in close agre ement with that of previous studies, e.g . E{B — V)= 0.88 


- 0.94 mag flMoffat et ah 


(jSuiatha fc Babu 


1979 ). 0.84 - 0.87 m ag (Chini fc Wink 


2006 ). and 0.76 - 1.00 mag f Jose et al.l 


19841) ■ 0.70 - 0.90 mag 


2011 ). The reddening of NGC 


1931 obtained from 14 early-type members ranges from E{B — V) = 0.51 - 1.01 mag. The 
mean reddening is {E{B — V)) = 0.74 ± 0.17 (s.d.) mag. This re sult is also commen surate 


with that of previous studies, e.g. E (B — V) = 0.49 - 0.93 mag (IV 


1.00 mag flPandev fc Mahra 


(iBonatto fc Bica 


I 


19861) . 0.55 - 1. 00 mag llBhat 


: et ah 


20091) . and 0.50 - 0.90 mag (iPandev et al.l 


offat et ah 


19791) . 0.47 


19941) . 0.52 - 0.72 mag 


2013al ). The dispersion in the 


reddening indicates that there is differential reddening across each cluster. The differential 
reddening in NGC 1624 appears to be less signihcant than that found in NGC 1931. 


Unlike early-type MS stars, it is difficult to determine the reddening for the PMS stars 


because the intrinsic colors can be altered by hot and cold spots on t 


activities, and the obscuration by a circumstellar disk fl Grankin et ah 


le sur face, accretion 


20071) . Simultaneous 


imaging (color) and spectroscopic (spectral type) observations are required to determine the 
reddening of individual PMS stars accurately. In the absence of such data, we estimated the 
reddening of PMS stars by using a weighted-mean reddening value at a given position from 
the reddening distribution of the early-type members, where the weight was given by an 
exponential function with respect to the distance from the individual early-type members. 
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The reddening law toward young open clusters provides an opportunity to study the 
spatial distribution of interstellar matter in the Galaxy as well as the dust evolution in 
SFRs. The ratio of total-to-selective extinction (Ry) is a basic tool to investigate the 
reddening law, being also a crucial parameter to correct the total extinction in the visual 
band. The parameter is closely related to the size distribution of dust grains. The Ry 


Greve 

2010; 

Hur et al. 

2OI2I 


20121) . Previous stuc 


reddening law {Ry = 3.1, 


Ghini fc Wink 


1984J: 


i es on NGC 1624 have adopted the normal 


Jose et al. 


2011 ). while 


Pandev et al. 


(l2013al) 


showed, from photometric and polarimetric data, that the reddening law toward NGG 1931 
is slightly different from that found in the general ISM. In order to check the reddening 
law toward these clusters we investigated vario us color excess ratios in the optical and IR 


pass 


j ands as performec 


2012 


Sung et al.l 


2013b 


in the series of stu dies (IKook et al. 


him et al. 


2014aUbh 


20101 : 


him et al. 


20111 


Hur et al. 


him et al 


fl2014aj) have noted that several mid - late-B-type MS stars in the young 


open cluster IG 1848 have a bluer color than that of normal B-type MS stars. Since the 
bluer color may not reflect the genuine photospheric color of the stars, the early-type 
members with U — B < 0.3 for NGG 1624 and with U — B < —0.2 for NGG 1931 were 
chosen to avoid such a color anomaly. We used 2MASS NIR data (tagged as a photometric 
quality flag ‘AAA’) as well as our optical data fo r the stars. For MIR data, the Spitzer 


GLIMPSE360 catalog flWhitnev et al 


2008 


201 ll) was used to obtain the R-MIR color 


excess. The various color excesses E(y — A) in the optical and IR passbands were computed 
by comparing the observed colors with the intrinsic relations from Paper 0 and Sung et 
al. (in preparation). The Ry is expressed a.s a linear relation for each color excess ratio 


E{V -\)/E{B-V) flGnetter fc Vrba 


1989 


Sung et al. 


2013bl ). We present the color excess 


ratios of the cluster members in Figure [9l The thick solid line (blue) in the hgure represents 
the slope corresponding to the normal reddening law {Ry = 3.1). 













































































21 


There is a foreground early-type star in the direction of NGC 1624 which was helpful 
to check the difference between the general ISM and the intracluster medium (ICM). As 
shown in the upper panels of Figure [HI all the color excess ratios consistently indicate 
that the ICM of NGC 1624 is indistinguishable from the general ISM. The mean Ry 
obtained from different color exce ss ratios was 3.12 ± 0.01 (s. d.). This result is in good 


agreement with that examined by 


Jose et ah 


(120111) , implying that dust evolution in the 


cluster had already progressed or no grain growth had occurred in the natal cloud. We 
adopted Ry = 3.1 and obtained a mean extinction of {Ay) = 2.85 ± 0.17 mag. On the other 
hand, there was no available foreground early-type star toward NGC 1931 for checking the 


foreground reddening law. Given that the reddeni ng law toward the neighbouring c 


NGC 1893 and S tock 8 turned out to be normal fjjose et ah 


2008 


Eswaraiah et ah 


usters 


2011 


him et ah 


2014bl) . we assumed that the foreground reddening law in the direction of NGC 


1931 was likely to be normal. In the lower panels of Figure [91 the color excess ratios of the 
members show a signihcant deviation from the normal reddening law. It implies that the 
size distribution of dust grain in the ICM of NGC 1931 is far different from that found the 
general ISM. Thus, the total extinction in V band toward the cluster should be divided into 
two components as below: 


Ay = fg X F(B - V)fg + Ry^ci X l^(B -V)- E{B - l/)fg] (2) 

where Ay, Ry,fg, Ry,ch B{B — G)fg, and E{B — V) represent the total V band extinction, the 
foreground Ry, intracluster Ry, foreground reddening, and total reddening, respectively. As 
mentioned above, Ry,{g = 3.1 was adopted for the foreground component. We simultaneously 
determined E{B — V)ig and Ry,c\ from the G-IR color excess ratios using a htting 
method. The star ID 800 shows signihcant excess emission at all the wavelengths. Another 
star ID 872 with a photometric quality hag ‘FEE’ in the 2MASS catalog also exhibits excess 
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emission in the MIR passbands. We did not include these stars in the examination of the 
reddening law. The mean foreground reddening was found to be E{B — I/)fg = 0.417 ±0.005 
(s. d.) mag, corresponding to a point of contact between foreground (thick solid line) 
and intracluster (thin solid line) components in the lower panels of Figure [H This value 


is in good agreement wit 
[E{B — V) = 0.418 mag, 


1 the sma 


Limetah 


lest re ddening of the neighbouring cluster NGC 1893 


2014b| . which is likely the amount of reddening in front 


of the clusters. The mean intracluster Rv,ci determined from R-IR color excess ratios was 
found to be 5.2 ± 0.3 (s. d.). The slope corresponding to the Rv,c\ is shown by thin solid 
lines in the lower panels of Figure O 

The E(y — I)/E{B — V) color excess ratio gives an intracluster Rv,c\ of 4.3 (dashed 
line). The value is significantly different from the mean Rv,c\ obtained from the R-IR colors 
excess ratios. We attempted to check our V — I colors. As discus sed in Section 2.3, th ere is 


no considerable systematic difference between the photometry of 


Pandev et ah 


f 2013a l and 


ours. The observations of NGC 1624 was made on the same night, and the multiwavelength 
study on the reddening law showed a consistent result in the optical - NIR passbands. 
Hence, our photometry and the adopted empirical relations of Paper 0 do not suggest any 
serious systematic problem. A similar aspect w as found in the reddening law toward the 


starburst cluster Westerlund 2 flHur et al.l 


201411 . The authors attributed this discrepancy 


to the unknown behaviour of the V — I color in heavily reddened situations. It needs to 
be confirmed whether or not the linear relation between Ry and E{V — I)/E{B — V) is 
applicable to extremely young and dusty SFRs. We leave this issue for future work in our 
survey project. 


Using the IR extinction curve of 


Fitzpatrick fc Massal (1200711 (equation 4) we tested 


the reddening law once more as shown in Figure dOl The extinction curve of NGC 1624 
(open circles in the upper panel) is close to the mean Galactic extinction curve (dashed line 
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^ ^iR = 1-06 and Ry = 3.00 - Fitzpatrick &: Massa 2007 ). while that of NGC 1931 still 
shows a conspicuous deviation from the mean curve with a somewhat large scatter. The Ry 
was determined by a fitting method, and the result was Ry = 3.8 ± 0.1 (1/A < l/rm“^). 
However, we cannot accept that a small specific column throughout the sky can be filled 
with significantly lar ge dust grains, although a global va riation of Ry along the Galactic 


longitude was found fjWhittet 


1977 : 


Sung fc Bessel 


2014) . We derived the IGM component 


in the extinction curve after subtracting a foreground reddening oi E{B — l/)fg = 0.42 
mag. The lower panel of Figure fTOl shows the extinction curve of the IGM component. 

We found Ry^c\ of 5.5 ± 0.3 using the same procedure as above. The result is reasonably 
consistent with that obtained from Figure [9] within the uncertainties. We note that the 
E{y — I)/E{B — V) color excess ratio seems to be ambiguous in this analysis because it is 
commensurate with the mean Galactic extinction curve as well as the IR extinction curve 
of the cluster {Ry^ci = 5.5). 

We adopted the segmented reddening law toward NGG 1931 as presented in Figure [9] 
[Ryfg = 3.1, E{B — 'F)fg = 0.42 mag, and Ry^ci = 5.2] and obtained a mean extinction of 
{Ay) = 2.97 ± 0.88 mag. This result is somewhat differe nt from that obtaine d through a 


polarimetric and photomet ric method {Ry = 3.2 - 3.3 


is that the Serkowski law fjSerkowski et al 


Pandev et ah 


2013al) . A reason 


19751 ) may not be applicable in the extremely 


young and dusty SFRs because of the efficiency and time scale for grain alignment, the 
probable complex structure of the magnetic held in the IGM, as well as the depolarization 
effect by foreground dust grains. The other reason is that the authors did not consider 
the foreground and intracluster components separately, and therefore their Ry could be 
underestimated. Nevertheless, their results also suggest that the reddening law in the 
direction of NGG 1931 is not normal. We conclude that the evolution of dust grains within 
the cluster is still in progress. 

















24 


3.5. Distance 


Young ope n clusters (< 10 Myr ) are a useful tracer of the local spiral arm structure 


in the Galaxy flDias &: Lepine 


20051 ). We have been using the zero-age main sequence 


(ZAMS) htting method to determine distances to open clusters in a series of studies. The 
canonical ZAMS fitting method should be made after correction for interstellar reddening. 
In order to reduce uncertainties arising from the reddening correction, Paper 0 introduced 
the reddening-independent indices as below: 


Qyj = V-2A5{V-I) 

(3) 

Qyj = V - 1.33{V - J) 

(4) 

Qyj, = V -1.17{V - H) 

(5) 

Qyj,^ = V-1.10{V-Ks) 

(6) 


The modified Johnson Q [Q' = (U — B) — 0.72(5 — V) — 0.0255(5 — VY] is used as the 
abscissa of new CMDs. These reddening-independent indices give a few advantages in 
ZAMS htting procedure. The indices are based on UBVIJHKs multicolor photometry, 
and thus the ZAMS htting to the CMDs of young open clusters can be made consistently 
with respect to diherent colors. The colors of early-type stars are less ahected by metallicity 
ehects because few strong metallic lines apart from a few light elements are, in general, 
found i n their spectra. W e have determined the distance of young open clusters using these 


indices (iLim et ah 


2014alJbh . 


A large fraction of high-mass stars constitute binary systems fjSana et al.l 120121 ). The 


luminosity of such stars appears higher than that expected from ZAMS stars at a given 
color or ehective temperature. In addition, physical properties, such as stellar rotation and 


overshooting, cause the MS band to be wider in the high-mass regime ( Martins fc Palacio; 
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2013 '). Therefor e, the faintest MS s 


as ZAMS stars fl.Tohnsnn fc Hiltner 


ars a t a given color in the CMDs have been regarded 


1956h . It is worth noting that some of the early-type 


members may possibly reveal a color anomaly as mentioned in the previons section. 
Allowing for this point we carefnlly fitted the ZAMS relations of Paper 0 to the lower ridge 
line of the clnster sequence in Figure ITT] after adjusting the relations above and below it. 
The lower ridge line could be conhned between the ZAMS relations (dashed lines) shifted 
by ± 0.2 - 0.3 mag for each cluster. The htted distance modulus of NGC 1624 and NGC 
1931 was 13.9 ± 0.2 mag (6.0 ± 0.6 kpc) and 11.8 ± 0.3 mag (2.3 ± 0.3 kpc), respectively. 
This result places NGG 1624 (Z = 155?356, b = 2?616) in the outer arm and NGG 1931 
{I = 173?898, b = 0?281) in the Perseus arm. 

Previous studies used either ZAMS or isochrone fitting in the determination of distance. 
The distance to NGG 1624 ob tained in this work is in good agreement with the results o f 


previous stud ies, e.g. 6.0 kpc 


and 6.1 - 6.5 flDeharveng et al. 


Moffat et al. 


1979 : 


Suiatha fc Bab 


20081) . It is interesting that 


2006: 


Jose et al. 


Suiatha fc Babu ( 2006 ) obtained 


201 ih 


a distance consistent with that of oth er studies although t 


systematic difference in V (Figure [2]). iGhini fc Winld (119841 ) obtained a distance of 10.3 kpc 


leir photometry shows a serious 


using spectroscopic parallax for 2 bright stars. In their study, the stars were classified as 
061 and BlIII, respectively. Since their spectral classification is so mewhat different from 


that of a recent study for the brightest star NGG 1624-2 (07f?p 


Sota et al. 


20 


been made in 


Walborn et al 


2010 


ll), their c 


Jose et al 


their previous results flPandev fc Mahra 


istance may be shifted systematic ally. A similar discuss ion has 


( 2011). In the case of NGG 1931 


distance of 2.3 ± 0.3 kpc. 


Bonatto fc Bica (2009) ob 


the 2MASS NIR photometry. Although 


986- 


Bhatt et al. 


Pandev et al. 


( 2013al ) refined 


19941) and published a new 


Moffat et al. 


ained a distance of 2.4 ± 0.3 kpc from 


(119791) reported a rather smaller value 


(1.8 kpc), our result is in good agreement with that of more recent studies. 
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4. HERTZSPRUNG-RUSSELL DIAGRAM 


The HRD is a basic tool to understand the evolution of stars. The effective temperature 
(Teff) of stars is a crucial parameter for constructing the HRD. We used the relations 
published in Paper 0 (table 5) to obtain the T^ff and bolometric correction for individual 
stars. The Tg// of the star NGC 1624-2 was determined from the spectral type- Teff 


relation. Although the 


05.5V (Moffat et ah 


Sota et al 


by 


201 


Wade et al 


classihcation fISota et al 


spec 


ral typ e and luminosity clas s of the s tar seem to be uncer tain 


1979h . 061 fIChini fc Wink 


198411 ■ 07f?p f Walborn et al. 


2010 


1, an d 06.5f?cp-08f?cp (a variation in the spectral type has been reported 


201211. we assumed the star to be an 07 MS star according to the recent 


2011 


Wade et al. 


201211 . It is worth noting that the absolute 


magnitude of NGC 1624-2 (My = -4.66 mag) is reasonably consistent with that of 07 MS 
stars (My = -4.90 mag) rather than that of 07 supergiant stars (My = -6.05 - -6.95 mag, 
see table 4 of Paper 0). The difference between an 07 MS and giant star is only about 
165 K in Teff, and the stars have the same bolometric correction value (Paper 0). If NGC 
1624-2 is a supergiant star, the difference between MS and supergiant stars increases up 
to 800 K in Tg// and 0.1 mag in bolometric correction. The T^ff of other early-type MS 
stars was inferred from the color- Teff relations. We averaged the Teff estimated from each 
color- '^eff relation with weights. 


We set the weights of Tg// estimated from U — B vs. T^fj relation to be 0.5, 1.0, 1.5, 
2.0, and 2.5 in the {U — B)q color range of 0.07 to -0.1 mag, -0.1 to -0.3 mag, -0.3 to -0.5 
mag, -0.5 to -0.7 mag, and -0.7 to -1.2 mag, respectively, set it to be 0 for late-type stars. 
The temperature sensitivity of the V — I color is not high enough to estimate the Teff of 
early-type stars. We assigned each weight of Teff inferred from B — V vs. Teff relation 
and V — I vs. '^eff relation to 1.0 and 0.5 for stars with {B — V)q < —0.24 mag. For MS 
stars with {B — V)o < 0.18 mag, the weights of Teff estimated from B — V vs. Teff relation 
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and V — I vs. T^ff relations were set to be 1.0 and 0.7, respectively. Since the V — I 


color is a goo c 


flBessell 


1995 


temperatnre indic ator for cool stars, we only nsed V — I vs. Tg// relations 


Bessell et al 


19981) for PMS stars. The bolometric correction valnes for all 


the members were inferred from their T^ff nsing table 5 in Paper 0. We present the HRD 
of NGC 1624 and NGC 1931 in Figure [121 


The brightest star in NGG 1624 is an 07 star flWalborn et al 


2010 


Sota et al. 


20111 ) 


If the star is a V 
to 


S star with solar metallicity, its lifetime is smaller than 10 Myr according 


MassevI (120131) . On the other h and, the spect r al type of the most luminous star in 


NGG 1931 is somewhat uncertain. 


Pandev et al. 


fl2013al) argued that the main ionizing 


sources of Sh2-237 may be two B2 stars. The MS lifetime of such stars is tens of millions 
of years. If the stars are in the MS stage, the lifetime is the upper age limit of the cluster. 
We superi mposed a few i sochro nes solid lines) constr ucted from the stellar evolution 


models of 


Ekstrom et al 


(l2012f) and 


Siess et al. ( 2000 ) on the HRD of NGG 1624 with 


several evolutionary tracks (dashed lines). The position of NGG 1624-2 in the HRD is well 


matched to the 4 Myr isochrone, and thus the turn-off age is a bout 4 Myr. We a 


the isochrones interpolated from the stellar evolution models of 


Brott et al. 


(12011 


so used 
I for the 


initial chemical composition of the Galaxy and Large Magellanic Glo ud (LMG), whe r e the 


evolutionary tracks with the similar initial rotation velocity to that of 


Ekstrom et al. 


( 120121 ) 


were selected. The isochrone for the chemical composition of the Galaxy gives an age of 3.3 
Myr, while the 3.8 Myr isochrone for the slightly lower LMG metallicity appears to well 
predict the Tg// and Mboi of NGG 1624-2. A systematic uncertainty of the turn-off age 
arising from the difference in metallicity may be about 0.5 Myr. 


In the case of NGG 1931, it is impossible to infer the turn-off age because the most 
luminous star ID 872 in NGG 1931 is below the ZAMS line. The 
star became bluer by ~ 0.1 mag from 2005 to 2013 (the data from 


U — B color of the 


Pandev et al. 


2013al) . 




















































28 


Another early-type stars ID 629 and 854 also showed a similar variation in the U — B color. 
The mysterions yonng stellar object Walker 90 fonnd in the yonng open clnster NGC 2264 
has also shows s uch a variation in its spectral type as well as colors for the last ~ 50 years 
(from A3 to B4, 


Perez et ah 


2008 


and references therein). The observational properties of 


the star, such as an abnormal reddening law of Ry = 3.6 ~ 7.0, the inverse P-Cygni profile 
of Balmer lines, IR excess emission, and non-photospheric UV conti nuum, indicate th at the 


star is likely an intermediate-mass PMS star with an accretion disk fjPerez et ah 


2008h . The 


star ID 872 in NGC 1931 also showed similar photometric properties (the variation inU — B 
color, an abnormal reddening law, and IR excess emission) to those of Walker 90 although 
its spectral features were not conhrmed. If the star is at the analogous evolutionary stage 


to Walker 90, the age of NGC 19 31 can be conjee 


within NGC 2264 (1.6 - 3.0 Myr, 


Sung fc Bessell 


ured from the age of the S Mon group 


2010h 


Most of the PMS members in the two clusters have masses smaller than 3 Mq. The 
number of PMS members identified in NGC 1624 is insufficient to investigate the age 
distribution. The star ID 1773, which is known as the main ionizing source of UCHII at 
the border of NGC 1624, is likely a PMS star at a very early evolutionary stage. The age 
of the star seems to be younger than 1 Myr, indicating that star formation is currently 


taking place in the cluster. We estimated the ma ss accretion rate o 


excess emission as performed in previous studies flPebnll et ah 


2000 


he star using its UV 


him et ah 


2014a 




The mass accretion rate of the young PMS star is M = 3.98 x 10~^Mq yr“^. On the other 
hand, more than half of the PMS members are likely to be younger than 4 Myr. There are 
three PMS members (2 Ha emission and 1 UV excess stars) near the ZAMS. These stars 
seem to be older than other members. However, if the stars have a nearly edge-on disk, 
their luminosities would be underestimated. In addition, accretion activities can lead the 
colors of the stars to be bluer, i.e. hotter Tg//. 




























-29- 

In the HRD of NGC 1931 (the right panel of Figure [12]), most of the PMS members are 
younger than 10 Myr, which was suggested by the open cluster data base WEBDA^. It is 
worth noting that WEBDA provided incorrect ages for some open clusters. The number of 
PMS members is relatively larger than that identihed in NGC 1624 because of its proximity 
as well as the variety of the available data. We investigated the age distribution of the 
PMS members (thick solid line) in Figure [T31 where the bin size is 0.5 Myr. In order to 
avoid binning effects, another histogram (thin solid line) was plotted by shifting the bins 
by half the bin size. A peak in the age distribution appears at 1.5 Myr. The median age is 
about 2.0 Myr with a spread of 4.5 Myr, where the spread was dehned as the age difference 
between the 10 and 90 percentiles in the cumulative age distribution of the PMS members 


fjSung fc Bessell 


20101) . Since all the PMS members were identihed using membership 


selection criteria based on the observational properties of warm circumstellar disks and 
accretion activities, the age distribution could be biased to young active stars. 


Previous studies e stimated the turn-off age of NGC 1624 to be 4 Myr fISuiatha &: Babu 


20061 : 


Jose et ah 


evolution models by 


201111. Our est i mate 


Brott et ah 


(4 Myr) is in good agreement with those. Other 


(120111) for Milky Way and LMC metallicities give a slightly 


younger age (3.3 and 3.8 Myr) . The upp e r valu e of the ages of PMS members (< 4 Myr) 


is also consistent with that of 


Jose et af 


f 201ll ). The turn-off age of NGC 1931 remains 


uncerta in because the evolutio nary stage of the brightest star is ambiguous. For this 


reason, 


Pandey et ah 


Bonatto &: Bical (12009!) obtained a somewhat larger turn-off age (10 ± 3 Myr), while 


(j2013al ) only suggested an upper value of 25 Myr. The latter reported a mean 


age of 2 ± 1 Myr for the PMS stars. This is in good agreement with the median age from 
this work. They also found a similar age spread (~ 5 Myr) from their CMD. We have 
further discussions on the age spread of PMS stars in Section 6. 


^ ht t p: / / WWW. univie .ac.at/webda/ 
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5. INITIAL MASS FUNCTION 

5.1. Derivation 


The membership selection of PMS stars is likely incomplete because the membership 
selection criteria used above are more likely to identify young active PMS stars rather than 
those in a quiescent phase. A statistical method was therefore used to hll in the incomplete 
mass bins of the IMF. 


We included all the stars within the PMS locus in the (V, V — I) CMD, that is, 
bona-hde members plus any held interlopers. These stars were then placed in the HRD 
following the same procedure as delineated in the previous section (see Figure [12]). The 
mass of the individual stars was estimated in the HRD by comparing their Tg// and Mi,oi 


wi th those from 
of Ekstrom et ah 


he evo lutionary tracks. For the MS members, the stellar evolution model 


(1201211 was used to estimate the mass. The mass of the PMS m embers and 


the held stars was inferred from the PMS evolution models of 


Siess et ah 


(I 2 OO 0 I) . In order to 


minimize the contribution of held interlopers, appropriate control helds were chosen within 
the observed regions (the shaded area in Figure [T|). The stars within the PMS locus in the 
control helds were assumed to be evenly distributed in the direction of each cluster. The 
apparent mass of these stars was estimated by using the same procedure as above, where 
the mean reddening and distance of the two clusters applied to the stars. We note that the 
apparent mass does not mean real stellar mass. 


The IMF is, in general, expressed as .^ = XkifiTs’ where A, Alogm, and S represent 
the number of stars within a given mass bin, the size of the logarithmic mass bin, and the 
area of the cluster, respectively. We applied a slightly large bin size of Alogm = 0.4 to 
include as many stars as possible within a given mass bin. An even larger bin size of 1.0 
was adopted for the highest mass bin [10 < M/Mq < 100 - we assumed the upper limit 
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of stellar mass (Mup) to be 100 Mq] because of the small number of high-mass stars. The 
number of stars was counted within each given logarithmic mass bin, and then that was 
normalized by the size of the mass bin and the area of the clusters. The apparent IMF of 
the held interlopers was obtained in the same way. In the case of NGC 1931, the extent 
of the control held was smaller than that of the cluster, and therefore the number of held 
interlopers was corrected by multiplying the areal ratio of the cluster to the control held. 
Finally, the IMF of the two clusters was obtained by subtracting the contribution of held 
interlopers. Figure [H] shows the IMF of the two clusters. In order to avoid binning ehects 
we shifted the mass bin by 0.2 and re-derived the IMF (open circle) in the same way. The 
error of the IMF was assumed to follow Poisson statistics. The arrows represent the lower 
limit of the IMF below the completeness limit of our photometry. 


Using a least-square method we determined the slope of the IMF in the mass range of 
3-27 Mq for NGC 1624 and 1-18 Mq for NGC 1931. The slope is F = —2.0 ± 0.2 for 


the former and F = —2.0 ±0.1 for the la; 


than the Salpeter/Kroupa IMF (jSalpeter 


er. 


1955 


le IMF of these clu sters appears far steeper 

1 as the results of 


Kroupa 


2001 


previous studies , e.g. F = —1 . 18 to —1.31 for NGC 1624 (jjose et ah 


2002h as we. 


for NGC 1931 f Pandev et ah 


201lh and P = -1.15 


2013al) . However, the small number of high-mass stars 


(> 10 Mq) and the upper limit of the stellar mass which was required for the normalization 


of the IMF may be responsible 
cluster NGC 2264 (dashed line. 


or the steep slope. T he IMF of the well-studied young open 


Sung &: Bessell 


2010l) is shown in Figure [H] for comparison. 


The IMF of the two clusters seems to resemble that of NGC 2264 within the uncertainties. 


According to a recent review on the IMF fjOffner et ah 


201411 the IMF of NGC 2264 is 


m 


good agreement with that of nearby SFRs. This may imply that the star formation in NGC 
1624 and NGC 1931 is not much different from that in the solar neighbourhood. More 
discussion on the slope of the IMF is addressed in the next section. 
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We estimated the total number of members and the masses of the two clusters by 
integrating their IMFs. The IMF in the low-mass regime (from the completeness limit down 
to 0.25 Mq) was assumed to be that of NGC 2264. NGC 1624 hosts 66llgJ members, and 
a total of 656^10^ stars constitute NGG 1931. The total mass of NGG 1624 and NGG 1931 
is about Mq and 510^7°^ Mq, respectively. The uncertainty in the total masses was 

propagated from the error of the IMF in given mass bins. The total mass is comparable to 


of af 


members from 


Sung & Bessell 

2010 

), 550+1° Mq 

for the Pleiades 

(Adams et ah 

2001) 


or Praesepe flKraus &: Hillenbrand 


of magnitude li ghter than t 


le most massive star 


(> 50,000 Mq, 


Glark et al.l 


2005 


Gennaro et ah 


jurst cluster Westerlu nd 1 in the Galaxy 


2011 


him et ah 


2013h . 


5.2. Implications of the Slope of the IMF 


The IMF of NGG 1624 and NGG 1931 shows a steep slope for stars with mass larger 
than 1 Mq. This could imply that low-mass star formatio n is dominant in the clusters. 


According to the metallicity variation in the Galactic disk flYong et ah 


20121) the chemical 


composition at the Galactocentric distance of NGG 1624 (i?GC = 13.7 kpc) and NGG 
1931 (i?GC = 10.3 kpc) is expected to be 0.1 - 0.3 dex lower than the solar metallicity. 
The temperature of a molecular cloud may therefore be higher in th e low metallicity 


envir onment because radiative cooling by metallic ions will be lower (iGaramazza et ah 


2 OI 2 I ). Furthermore, there is a rela tion between the mass accretion rate of a star and the 


temperature of a cloud (M oc 


Susa et ah 


20141) . As a result, high-mass stars may be 


formed advantageously in such a low metallicity environment, rather than active low-mass 
star formation. However, the IMF of the two clusters is inconsistent with this expectation. 


In order to properly interpret the IMF of the clusters we will hrst discuss a few issues 
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concerning the snbtraction of field stars, dynamical evolution and stochastic effects. 


We carried out checks to see whether the subtraction of held interpolators had been 
adequately carried out because had the contribution of held interlopers been improperly 
subtracted within a given mass bin, the resultant IMF would be misleading. We examined 


the V band luminosity function (L 
population model for the Galaxy 


) of stars in each control held using the synthetic stellar 


Robin et al. 


20031) . The synthetic LF was scaled to 


have the same area as that of each control held. The observed LF was compared with the 
synthetic one. The LF of stars in the control held for NGC 1931 (/ = 173?90, b = 0?28) is in 
good agreement with the synthetic one, while the stellar population model overestimated 


the number of stars toward 


at I = 90° (IRobin et ah 


VGG 1624 at R ~ 19 mag. Because of the Galactic warp found 


200311 ■ the synthetic stellar population in the direction of NGG 1624 


{I = 155?356, b = 2?616) which is linearly interpolated from those between I = 90° and 180° 
may not warrant consistency with the observational result. Alternatively, we obtained the 
IMF of NGG 1624 selecting diherent control helds (the west region Aa = —6f0, A5 = +3f5 
or the south-west region Aa = —7'5, A6 = —6f5 in Figure [I]). The slope of the IMF was 
about F = —2.2 ± 0.3 and —2.0 ± 0.3, respectively. Hence, the uncertainty arising from the 
subtraction of held interlopers has a negligible contribution to the slope. 


Kroupal (120011) has presented a segmented IMF based on compiled data. In his 


alpha-plot, a non-negligible scatter was found in the mass range of 1 - 100 Mq. The 
author attributed the cause of the scatter to dynamical ehects and Poisson noise. It used 
to be believed that young open clusters were too young to have undergone dynamical mass 
segregation. Howeve r, a few young ope n clusters sho w clear signs of m ass segregation. 


e.g. Arches cluster (IHabibi et ah 


(Raboud fc Mermilliod 


Westerlund 1 flGennaro et al 


1998 


2013h. NGG 3603 flPang et al. 


Sung et al. 


2011 


2013 


Limetah 


2013h. NGG 6231 


, ONG (IHillenbrand fc Hartmann 


1998h . 


20131) . etc. Several studies, through multiple 
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sets of numerical simul ation, suggested t 


has a dynamical origin flMcMillan et al. 


that the mass segregation 

bund in voung clusters 

2007: 

Moeckel & BonnelJ 

2009; 

Allison et al. 

2009 


2010h . If dynamical evolution had efficiently operated within young clusters, a number 


of low-mass members would have evaporated into the outskirts of the cluster region. 
Consequently, one would expect the IMF to be biased toward bright stars, and therefore 
show a shallow slope. The steep slope of the IMF in this work is inconsistent with this 
expectation of the effects. Thus, the contribution of the dynamical effect may be less 
signihcant for the clusters. 

The total number of members belonging to NGC 1624 and NGC 1931 is so small, that 
its measured IMF is vulnerable to stochastic effects. In order to test for this effect on the 
slope of the IMF, multiple sets of Monte Carlo simulation were conducted for three cases 


with the tota 


number of cluster members (iV* = 100, 1,000, and 3,000). The Kroupa IMF 


flKroupall200ll) was used as the underlying IMF for stars with mass larger than 0.5 Mq. 


The simulation for each case was made 10,000 times. The IMF of the model clusters was 
derived using the same procedure described as above. The size of the mass bin was set to 
A logm = 0.4, while a bin size of 1.0 was adopted for the highest mass bin if the number of 
stars with mass larger than 10 Mq was smaller than 5. Figure [15] displays the distribution 
of F with the total number of cluster members. The dispersion of the distribution appears 
to be large in model clusters with a small number of members, while the IMF of model 
clusters with a large number of members shows a higher consistency with the underlying 
IMF. However, the possibility that the IMF of a given cluster has a slope steeper than 
F = —2.0 is only about 10 percent in the case of iV* = 100, and the likelihood decreases 
with the increase of the total number of cluster members. It implies that the possibility 
that the IMF of the two clusters in the present work exhibits a steeper slope than —2.0 is 
about 1 percent. Therefore, the stochastic effect seems not to explain our results reasonably. 
Although the Kroupa/Salpeter IMF provides a representative slope of the stellar IMF in a 
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given mass range, it remains to be seen whether the slope is universal or not. Future work 
on the IMF of small, young open clusters will observationally be used to examine this issue. 


DISCUSSIONS ON THE AGE SPREAD OF PMS STARS 


Falla fc Stahlerl fll999l. l2000l. |2002[) argued that the star formation activity within 
nearby SFRs has persisted for about 10 Myr and that the star form ation rate rnay be 


accelerating in the present epoch, on timescales of 1 - 2 Myr. Later, 


Falla et al. 


fcnobh 


found a signihcant depletion of lithium in four FMS members of the Orion Nebula Cluster 
(ONC), and the result supported their idea because it is known that most of lithium at 
the surface of stars with subsolar m ass are destroyed on timescales of ~ 10 Myr. On the 


other hand. 


Hartmannl (12001 


20031) refuted the extended star formation timescale invoking 


the influences of binarity, variability, extinction, accretion activ ities, the inclusion o f held 


interlopers, and the birthline effects of intermediate-mass stars. 


Jeffries et al. 


( 20071) found 


a small age spread of 2.5 Myr for NGC 2169. There are still many debates on the age 
of FMS stars although the age distribution is thought to be a key to understanding the 
processes of cluster formation. 

Only a small number of FMS stars were identihed as members of NGC 1624, so we 
were unable to statistically investigate the age spread. The number of the FMS members 
found in NGC 1931 barely allowed us to study their age spread, although the membership is 
still incomplete. The age distribution of the FMS members shows a spread of 4.5 Myr. We 
attempted to interpret the observed age distribution. Several observational uncertainties 
affecting the Te// and luminosity of PMS stars were co nsidered, e.g. photometric errors. 


reddening correction, and variability. 


Hartmannl (120031) pointed out that the inclusion of 


held interlopers can lead the observed age spread to be larger. However, we believe that 
most of the PMS stars identihed in this work are bona-hde members of NGC 1931 because 
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the membership selection criteria were based on the well-known observational properties 
of such stars together with the fundamental parameters of the cluster (see Section 3.2). 
Contamination by held interlopers was excluded as a consideration. 


We assessed the contribution of the observational uncertainties to the age spread 
of PMS stars using a few sets of simple Monte-Carlo simulations. NGC 1931 probably 
comprises more than 600 stars according to the integration of the IMF. We set the number 
of artihcial PMS stars to 600. The underlying IMF was assumed to be our result (F = —2.0) 
in the mass range of 0.5 to 5 Mq. The age of the stars was set to be 1.5 Myr (the peak 
age of NGC 1931 ) with no intrinsic age spread. Based on the PMS evolutionary models of 


Siess et ah 


(I 2 OOOI) artihcial PMS stars were generated. T^ft and lu r ninosi ty we re converted 


i nto V n magnitude and (V — /)o color using the relations of iBesselll (1l995l) and 


Bessell et ah 


fjl998h . We introduced photometric errors with a similar distribution to the observed one. 
The distribution of reddening was assumed to be a normal distribution, where the mean 
reddening and 3cr dispersion were set to be 0.74 and 0.17 mag, respectively. In order to 
reproduce the extinction in the V band we adopted the segmented reddening law as shown 
in Section 3.4 [E{B — V){g = 0.42 mag, Rv,fg = 3.1, and Rv,ci = 5.2]. We re-analyzed the 
photometric data of the artihcial stars to investigate their age spread after correcting for 
the mean reddening. Figure [I6] displays the result of the Monte-Carlo simulation. The 
resulting age spread of the artihcial stars was about 0.8 Myr although the age of a few 
PMS stars was overestimated. Hence, the photometric errors and the uncertainty in the 
reddening correction may not be enough to explain the observed age spread (4.5 Myr). 


The color variability of PMS stars can also increase the age spread because the 
position of the stars in the HRD is very sensitive to Te//- We investigated the variation 
of their physical quantities in the HRD using the data observed at two diherent epochs. 
The variable PMS stars identihed in this work changed by -1.2 - 1.5 mag in V and -0.5 
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0.5 mag m V — I over the 1.2 years. The variation in V appeared larger than that 


related to accret 

ion ac 

divities, as we 

. as to 

Basri & Batalha 

1990; 

Grankin et ah 

2007) 


as to variable reddening 


20071). Figure [m shows 


the HRD of the variable PMS stars. Assuming that the variation of reddening in the 
line of sight is negligible on a timescale of 1 year, we examined the change in age. The 
intermediate-mass PMS star ID 727 (the most luminous star in the HRD) varied the most 
in Tgff and luminosity over the 1.2 years. As a result of this variation, the apparent age 
changed by 2 - 3 Myr. Some low-mass stars also revealed as large an age variation although 
the change in the physical quantities was not as large as those of the star ID 727. This is 
because the resolution for distinguishing the age of PMS stars becomes low in the low-mass 
regime of the HRD. Consequently, the variability of PMS stars very likely contributes to 
the observed age spread. 


The fraction of PMS variables is at least 4.6 percent (16/351), where the total number 
of PMS stars (denominator) down to our detection limit {V ~ 22 mag) was estimated from 
integration of the IMF. This value is likely a lower limit. We expect more PMS variables to 
exist in the cluster because stars varying on much shorter or longer timescales may not have 
been identihed in our observations. A couple of Monte-Carlo simulations were performed 
assuming that 30 and 50 percent of the PMS members are variable. The variation in V and 
V — I was set to have a normal distribution with the same dispersion as the observed one. 
The photometric errors and reddening were treated with the same procedure as above. If 
we assume a variable fraction of 30 percent, the resultant age spread was about 1.5 Myr. A 
spread of 2.6 Myr was found for a variable fraction of 50 percent. These simple simulations 
showed that the variability of PMS stars can contribute a portion of the observed age 
spread. However, additional sources may still be required to explain the observed spread of 
4.5 Myr. 
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We briefly mentioned in Section 3.3 that the structure of NGC 1931 may have resulted 
from its star formation history. If star formation propagated in a specihc direction, one 
could expect there to be a correlation between the location and age of PMS stars. A 
systematic age variation on a few parsec scale has been found in the yo ung open cluster 


NGC 1893 I Sharma et ah 


29971: 


Pandey et af 


2913b 


Tnm et ah 


2914m . Such a star 


formation history can influence the age spread of PMS stars. We investigated the age 
variation along the declination axis, as the cluster is elongated in the north-south direction 
(Figure [8]). PMS members with masses larger than 1.5 Mq were excluded from this analysis 
because evolutio nary models for intermediate-mass PMS stars are likely to overestimate the 


age of the stars fjSung et ah 


19971. 


2994J: 


Hartmann 


19991). We considered only stars younger 


than 4 Myr to find their age variation clearly. Figure [18] shows the age variation of the PMS 
members with declination. The age appears to decline with the increase in declination. It 
implies that star formation has propagated from south to north with a projected velocity 
of 6.8 km s“^ according to the slope of the age variation. The maximum age difference 
between the stars in the northern and southern groups is about 9.7 Myr, on average. Hence, 
the star formation history in NGC 1931 seems to contribute a small portion of the observed 
age spread. In conclusion, a genuine age spread of PMS members is required to explain the 
observed age spread, however the size of the sprea d may not be as large as the timescale 


19 Myr) argued by 


Palla fc Stahlerl (2999 


299211 


7. SUMMARY 

NGC 1624 and NGC 1931 are small, young, open clusters in the direction of the 
Galactic anticenter. These clusters are helpful in studying the property of star formation 
in a different environment from the solar neighbourhood. We performed UBVI and Ha 
photometry for the clusters as part of the SOS project. This work provided homogeneous 
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photometric data as well as the comprehensive results for the clusters. 

The members of the clusters were selected using various criteria. For NGC 1624, we 
found 14 early-type and 14 PMS members (3 variable, 3 UV excess emission, and 9 Ha 
emission stars and candidates) from photometric criteria and Ha photometry. On the other 
hand, 14 early-type and 71 PMS stars (16 variable, 4 UV excess emission, 54 MIR excess 
emission, and 27 Ha emission stars and candidates) were identihed as the members of NGC 
1931 using photometric criteria, Ha photometry, and the GLIMPSE360 catalog. 

The reddening of individual early-type members was determined from the {U—B, B — V) 
TCD. The mean reddening was {E{B — V)) = 0.92T0.05 mag for NGC 1624 and 0.74±0.17 
mag for NGC 1931. The reddening law toward the clusters was examined using various 
color excess ratios from the optical to NIR passbands (to MIR passbands for NGC 1931). 
We conhrmed that the reddening law toward and in NGC 1624 is normal. On the other 
hand, the early-type members of NGC 1931 exhibited color excess ratios that far deviated 
from the normal reddening. It implies that the size distribution of dust grains in the NGC 
1931 ICM is somewhat different from that found in the general diffuse ISM. The ratio of 
total-to-selective extinction for the ICM of NGC 1931 was estimated to be Rv,c\ = 5.2 ± 0.3, 
indicating that the evolution of dust grains is still ongoing. 


We carried out ZAMS htting to the reddening-independent CMDs of the early-type 
members. As a result, NGC 1624 was found to be 6.0 ± 0.6 kpc far away from the Sun, and 
probably located in the outer arm. The distance to NGC 1931 was about 2.3 ± 0.2 kpc, 
which places the cluster in the Perseus arm. Reddening-corrected CMDs wer e converted 


Bessell et ah 

1998; 

Sune et ah 

2013a1 


temper ature scales and bolometric correction (IBessell 


1995 


2013aJ). We performed isochrone fitting to the HRD of the 


clusters. A turn-off age of NGC 1624 was estimated to be 4 Myr, and the PMS members are 
likely to be younger than 4 Myr. The star ID 1773 associated with the UCHH region turned 
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out to be a very young PMS star with a high mass accretion rate of 3.98 x 10 “®Mq yr“^. 
The most luminous star in NGC 1931 seems to be a very young star with an accretion disk 
given that the photometric properties were similar to those of Walker 90 in NGC 2264. 
Thus, the age of the star is likely to be 1.6 - 3.0 Myr. We found a median age of 2 Myr 
with a spread of 4.5 Myr for the age distribution of the PMS members. 


Finally, we derived the IMF of the clusters and determined the slope in the high-mass 
regime (1 > Mq). The slope was F = —2.0 ± 0.2 for NGC 1624 and F = —2.0 ± 0.1 for 
NGC 1931, respectively. These results were somewhat different from the Salpeter/Kroupa 
IMF. We attributed the steeper slope to a stochastic effect, the small number of stars in the 
highest mass bin (> 10 Mq), and the assumed upper limit of stellar mass. The IMF of the 
cluste rs exhibits a very sim ilar shape to that of the well-studied young open cluster NGC 


2264 fjSnng fc Bessell Il2010l ). It implies that the property of star formation in the clusters 


may not be far different from that in nearby SFRs. 
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Fig. 1.— Finder charts of NGC 1624 (left) and NGC 1931 (right). Stars brighter than G = 18 
mag are plotted, and the size of the circles is proportional to the brightness of individual stars. 
The positions of stars are relative to the 0-type star NGG 1624-2 {a = 04^ 40™ 37'^.3, 6 = 
+50° 27' 41".l) for NGG 1624 and to BD+34 1074 (a = OS'* 31™ 26h4, 6 = +34° 14' 43".0) 
for NGG 1931, respectively. Squares outlined by blue solid lines represent the region observed 
by the Mont4k GGD camera, and the other square (green dashed line) shows the held of 
view of SNUGam. The shaded regions outline the control helds, which are used to estimate 
the density of held interlopers. 
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Fig. 2.— Comparison of our photometry with a few previous sets of CCD photometry. The 
upper and lower panels represent the difference between photometric datasets with respect 
to the V magnitude, and the V — I, B — V, and U — B colors for NGC 1624 (1) and NGC 
1931 (2). Bold dots (black) displays the consistency between our data obtained with two 
different observatio nal setups. Squa r es (b lu e), triangle s (red) , and crosses (green ) denot e 


the photometry of 


Suiatha &: Babu ( 2006 1. 


Jose et ah 


(1201111 . and 


Pandev et af 


(l2013all . 


respectively. Stars brighter than V = 17 mag are used in the comparison to avoid the large 
photometric errors of faint stars. 
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Fig. 3.— The radial surface density profile of NGC 1624. The error at a given surface density 
is assumed to follow Poisson statistics. Solid and dashed lines represent the mean surface 
density of field stars and its standard deviation, respectively. The surface density continues 
to decrease until it reaches the mean surface density of field stars. The radius of the cluster 
is about 2^5 (arrow). 
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Fig. 4.— Color-color diagrams of NGC 1624 (upper) and NGC 1931 (lower). The small dots 
represent all the observed stars. Different symbols identify early-type main sequence mem¬ 
bers (black bold dots), Ha emission stars (red circles), variables (blue triangles), ultraviolet 
excess stars (violet pluses). X-ray emission objects (large blue crosses), an X-ray emission 
candidate (a small blue cross), and mid-infrared excess stars (green squares). Since most 
pre-main sequence stars are too faint to be detected in the U band, they do not appear in 
the {U — B, B — V) diagram. The intrinsic and reddened color-color relations are overplotted 
by a solid line (blue) and dashed line (red) in the left-hand side panels. The mean reddening 
of {E{B — V)) = 0.92 mag for NGG 1624 and 0.74 mag for NGG 1931 is adopted for the 
reddened relations. The dotted line (cyan) denotes a hducial line to identify UV excess stars. 
In the right-hand side panels, the solid line represents the empirical photospheric level of 
normal main sequence stars, while the dashed and dotted lines are the lower limit of Ha 
emission stars and Ha emission candidates. 
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Fig. 5.— Color-magnitude diagrams for NGC 1624 (upper) and NGC 1931 (lower). Left- 
hand side panels: Id —/ vs. V diagram. Dashed lines (cyan) represent the empirical pre-main 
sequence locus. Middle panels: B — V vs. V diagram. Right-hand side panel: U — B vs. V 
diagram. The solid lines are the reddened zero-age main sequence relations of Paper 0. The 
arrow displays a reddening vector corresponding to Ay = 3 mag. The other symbols are the 
same as in Figure HI 
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Fig. 6.— The ([3.6] — [4.5], V — I) color-color diagram for NGC 1931. Dashed line represents 
the lower limit of the mid-infrared excess stars. The other symbols are the same as in 
Figure 01 



V-I 





54 


0.5 

•V 

0.4 

0.3 

0.2 


NGC 1624 


0 . 


0.0 


O © O O ©•©■©■©" 


0 


2 


14 16 18 

V magnitude 


20 



0 . 


Fig. 7.— Variability test in the V band. The upper and lower panels show the photometric 
errors of stars observed in NGC 1624 (upper) and NGC 1931 (lower). The red circles and 
vertical bars indicate the mean and standard deviation of the photometric errors for a given 
magnitude bin. Triangles denotes variable candidates. 
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Fig. 8.— The spatial distribution of members of NGC 1624 (left) and NGC 1931 (right). 
The solid, dashed, and dot-dashed lines show contours corresponding to 80, 60, and 40 
percent level of the maximum surface density. The size of each symbol is proportional to the 
brightness of individual stars. The other symbols are the same as in Figure 01 
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Fig. 9.— Color excess ratios obtained from the early-type star toward NGC 1624 (upper 
row) and NGC 1931 (middle and lower panels). The dots and the open circle denote the 
cluster members and a foreground early-type star. Thick solid (blue), thin solid (green), and 
dashed (cyan) lines corresponds to the normal reddening law {Ry = 3.1) and the intracluster 
reddening relations (i?yci = 5.2 and 4.3). The extinction toward NGC 1624 follows the 
normal reddening law, while the color excess ratios toward NGC 1931 exhibit a deviation 
from the normal reddening law. See the main text for details. 
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Fig. 10.— Infrared extinction cnrves of NGC 1624 and NGC 1931. Upper panel; color excess 
ratios in the direction of NGG 1624 (open circle) and NGG 1931 (bold dot). Lower panel: 
color excess ratios for the intraclu ster medinm of NGG 1931. T he dashed lines represent the 


mean Galactic extinction curve of 


Fitzpatrick &: Massal 1 20071) {Ry = 3.00 and km = 1.06). 


Fitting solutions are drawn by solid lines (upper - Ry = 3.8 and km = 1.5, lower - Ry^ci = 5.5 
and km = 2.6). 








Fig. 11.— Zero-age main sequence fitting to the early-type MS me mbers of NGC 1624 


(upper) and NGC 1931 (lower). The zero-age main sequence relations of 


Sung et ah 


ll2ni3ah 


were htted to the lower ridge line of the members. The solid lines (blue) represent the 
adopted distance moduli of 13.9 ± 0.2 (6.0 ± 0.6 kpc) and 11.8 ± 0.3 mag (2.3 ± 0.3 kpc) for 
NGC 1624 and NGC 1931, respectively. The dashed lines are the ZAMS relations adjusted 
by the htting errors. 
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Fig. 12.— The Hertzsprung-Russell diagram of NGC 1624 (left) and NGC 1931 (right). A 
few isochrones (solid line) with different age (0., 1, 1.5, 4, a nd 10 Myr) are superimposed on 


the d iagram with several evolutionary tracks (dashed line, 


Siess et af 


2000 


Ekstrom et ah 


2 OI 2 I ). The other symbols are the same as in Figure 01 
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Fig. 13.— The age distribution of PMS members identified in NGC 1931. The histograms 
are based on different binning of the same stars. A peak age appears at 1.5 Myr, and the 
median age is about 2 Myr with a spread of 4.5 Myr. 
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Fig. 14.— The initial mass function of NGC 1624 (left) and NGC 1931 (right). To avoid 
the binning effect we shifted the mass bin by 0 .2 and re-derived the I MF (open circle) using 


the same procedure. The IMF of NGG 2264 ( Sung fc Bessell 


i]) 


201011 is shown as a dashed 


line with its uncertainty (shaded area) for comparison. Arrows denote the IMF below the 
completeness limit of our photometry. See the main text for details. 
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Fig. 15.— Distribution of the IMF slope (F) with the total number of cluster members. A 
Monte-Carlo Simulation was performed 10,000 times for the model clusters consisting of 100, 
1,000, and 3,000 members, respectively. The histograms are based on different binning of 
the same artihcial stars (thick and thin solid lines). 
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Fig. 16.— The color-magnitude diagram (left), Hertzsprung-Russell diagram (middle), and 
the age distribution (right) of artihcial stars. This simulation was conducted to reproduce 
the age distribution of PMS members in NGC 1931. The solid line and arrow in the left-hand 
panel denote the empirical ZAMS relation (Paper 0) and the reddening vector, respectively. 
In the middle panel, the solid and dashed lines represent isochrones (0 and 1.5 Myr) and 


evolutionary tracks for PMS stars fjSiess et ah 


2000 


Ekstrom et ah 


20121) . The histograms 


in the right-hand panel are based on different binning of the same stars (thick and thin solid 
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Fig. 17.— Variable stars in the Hertzsprung-Russell diagram of NGC 1931. Open and filled 
circles represent the position of the variable PMS stars observed on 2011 October 29, and 
on 2013 February 5, respectively. Straight solid lines trace the variations. Several PMS 


evolutionar y tracks and iso c 


solid lines fjSiess et ah 


2000 


irones for 0, 1, 3, 5 , and 10 Myr are indicated by dashed and 


Ekstrdm et al. 


201211 
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Fig. 18.— The age variation of the PMS members of NGC 1931 along the declination axis. 
The mean age and standard deviation at a given declination bin are shown as an open circle 
and error bar. The dashed line indicates the peak age of 1.5 Myr in the age distribution 
(Figure [13]). 
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Table 1. Observation Log 


Target 

Date 

Filter 

Exposure time 

Maidanak Astronomical Observatory 




NGC 1624 

2006.11.24. 

I 

3s, 60s 



V 

3s, 180s 



B 

5s, 300s 



U 

15s, 600s 

Steward Observatory 




NGC 1624 

2011.10.29. 

I 

5s, 60s 



V 

5s, 180s 



B 

7s, 300s 



U 

15s, 600s 



Ha 

30s, 600s 

NGC 1931 

2011.10.29. 

I 

5s, 60s 



V 

5s, 180s 



B 

7s, 300s 



U 

15s, 600s 



Ha 

30s, 600s 


2013.02.05. 

/ 

5s, 180s 



V 

5s, 180s 



B 

7s, 300s 



U 

30s, 600s 
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Table 2. 


Atmospheric 


Extinction Coefficients and Photometric 


Zero points 


Date 

filter 

ki 

^2 

C (mag) 

2006.11.24. 

I 

0.056 ± 0.013 


23.229 ± 0.008 


V 

0.118 ±0.014 


23.632 ± 0.010 


B 

0.231 ± 0.003 

0.023 ± 0.002 

23.498 ± 0.010 


U 

0.403 ± 0.002 

0.018 ±0.001 

21.727 ±0.008 

2011.10.29. 

I 

0.046 ± 0.006 


22.196 ± 0.010 


V 

0.117 ±0.006 


23.543 ± 0.012 


B 

0.250 ± 0.005 

0.026 ± 0.002 

23.572 ± 0.009 


U 

0.460 ± 0.007 

0.020 ± 0.003 

22.109 ± 0.011 


Ha 

0.094 ± 0.009 


19.549 

2013.02.05 

I 

0.045 ± 0.008 


22.170 ± 0.009 


V 

0.120 ±0.008 


23.560 ± 0.007 


B 

0.232 ± 0.008 

0.023 ± 0.002 

23.548 ± 0.006 


U 

0.444 ± 0.018 

0.031 ±0.005 

22.069 ± 0.008 








- 68 - 


Table 3. Comparison of Photometry 


Paper 

AV^ 

N(m)^ 

A(V - I)^ 

N(m)^ 

A(S-V)1 

N(m)^ 

A([/-S)l 

N(m)^ 

NGC 1624 









This (MAO - SOoct) 

0.003 ± 0.013 

121(5) 

0.007 ± 0.018 

124(2) 

-0.006 ± 0.023 

123(3) 

0.001 ± 0.035 

120(5) 

Suiatha Babu f2006'l 

-0.872 ± 0.018 

29(5) 

-0.046 ± 0.014 

29(5) 

-0.032 ± 0.022 

31(3) 

-0.073 ± 0.064 

31(3) 

Jose et al. (^20^^^ 

0.020 ± 0.034 

125(5) 

0.031 ± 0.026 

180(10) 

0.017 ± 0.041 

129(1) 

-0.063 ± 0.077 

106(13) 

NGC 1931 









This (SOpeb - SOoct) 

0.004 ± 0.011 

92(15) 

0.013 ± 0.017 

102(5) 

-0.011 ± 0.018 

96(11) 

0.004 ± 0.027 

97(7) 

Pandev et al. f2013a) 

-0.023 ± 0.036 

142(12) 

-0.012 ± 0.022 

136(20) 

0.014 ± 0.022 

141(13) 

0.023 ± 0.063 

136(1) 


^This - Others 


N and m indicates the number of compared stars and excluded stars 










Table 4. Photometric Data for NGC 1624 


ID 

«^J2000 

^J2000 

+ 

I 

V - I 

B -V 

U - B 

H - C “ 

ey 

ej 

eV-J 

^B-V 

^U-B 

^H-C 


^obs 



2MASSID 

Ho 

2325 

04 

40 

39.03 

+50 

34 

39.2 

21.479 

19.533 

1.931 

1.367 



0.091 

0.086 

0.125 

0.194 



1 

1 

1 

1 

0 

0 

. 

. 

2326 

04 

40 

39.05 

+50 

28 

24.4 

18.969 

17.068 

1.883 

1.625 

1.419 

-0.353 

0.010 

0.002 

0.023 

0.056 

0.176 

0.057 

3 

5 

3 

2 

1 

2 

04403905+5028245 

h 

2327 

04 

40 

39.07 

+50 

27 

15.1 


18.484 






0.030 





0 

2 

0 

0 

0 

0 

- 

- 

2328 

04 

40 

39.08 

+50 

27 

13.4 


19.119 






0.226 





0 

2 

0 

0 

0 

0 

- 

- 

2329 

04 

40 

39.08 

+50 

30 

41.7 

21.374 

19.444 

1.858 

1.138 



0.085 

0.004 

0.060 

0.113 



2 

3 

2 

1 

0 

0 

- 

- 

2330 

04 

40 

39.10 

+50 

25 

35.9 

19.491 

17.965 

1.518 

1.261 

0.658 

-0.145 

0.001 

0.001 

0.005 

0.018 

0.054 

0.061 

2 

3 

2 

2 

2 

1 

04403910+5025356 

- 

2331 

04 

40 

39.12 

+50 

27 

27.6 

20.972 

18.660 

2.256 

1.732 



0.006 

0.001 

0.040 

0.162 



2 

3 

2 

2 

0 

0 

04403915+5027275 

- 

2332 

04 

40 

39.19 

+50 

27 

05.9 


19.340 






0.042 





0 

2 

0 

0 

0 

0 

- 

- 

2333 

04 

40 

39.20 

+50 

27 

07.9 

21.751 

18.229 

3.523 




0.186 

0.012 

0.186 




1 

2 

1 

0 

0 

0 

04403919+5027075 

- 

2334 

04 

40 

39.22 

+50 

27 

40.2 

16.956 

15.828 

1.138 

0.844 

0.251 

-0.327 

0.003 

0.012 

0.016 

0.010 

0.011 

0.007 

4 

5 

4 

4 

3 

2 

04403923+5027402 

H 

2335 

04 

40 

39.24 

+50 

18 

01.8 

20.933 

18.638 

2.294 

1.798 



0.083 

0.007 

0.083 

0.207 



1 

2 

1 

1 

0 

0 

04403923+5018016 

- 

2336 

04 

40 

39.24 

+50 

27 

14.2 


18.519 






0.027 





0 

2 

0 

0 

0 

0 

- 

- 

2337 

04 

40 

39.26 

+50 

19 

06.5 

21.818 

20.007 

1.811 




0.152 

0.145 

0.210 




1 

1 

1 

0 

0 

0 

- 

- 

2338 

04 

40 

39.30 

+50 

27 

11.6 

15.843 

14.918 

0.920 

0.691 

-0.039 

0.219 

0.001 

0.003 

0.003 

0.005 

0.005 

0.006 

4 

5 

4 

4 

4 

2 

04403929+5027li6 

- 

2339 

04 

40 

39.32 

+50 

28 

47.4 

18.428 

17.077 

1.350 

1.059 

0.428 

0.156 

0.013 

0.023 

0.011 

0.008 

0.020 

0.041 

4 

5 

4 

3 

2 

1 

04403931+50284© 

- 

2340 

04 

40 

39.32 

+50 

31 

39.1 

21.388 

18.611 

2.801 

1.465 



0.063 

0.006 

0.045 

0.244 



2 

3 

2 

1 

0 

0 

04403932+50313^2 

- 

2341 

04 

40 

39.34 

+50 

27 

19.6 

13.204 

12.338 

0.859 

0.623 

-0.298 

0.147 

0.005 

0.008 

0.001 

0.001 

0.001 

0.018 

3 

5 

3 

3 

3 

1 

04403935+5027196 

- 

2342 

04 

40 

39.38 

+50 

33 

01.5 

16.611 

15.251 

1.352 

1.131 

0.627 


0.004 

0.002 

0.004 

0.004 

0.014 


2 

3 

2 

2 

1 

0 

04403938+5033015 

- 


— C represents the Hex index [= Hex — (V + I)/2] 


H: Hex emission stars; h: Hex emission star candidates 







Table 5. Photometric 


ID 

Q:J2000 

^J2000 

y 

I 

V - I 

B - V 

U - B 

H - C 

705 

05 

31 

23.46 

+34 

12 

27.6 

21.091 

18.103 

3.031 



-1.609 

706 

05 

31 

23.47 

+34 

10 

47.2 

21.032 

18.326 

2.690 

1.529 


-0.555 

707 

05 

31 

23.47 

+34 

13 

28.0 

19.460 

17.885 

1.567 

1.205 

0.593 

0.038 

708 

05 

31 

23.48 

+34 

12 

18.4 

21.762 

18.790 

2.950 




709 

05 

31 

23.49 

+34 

13 

06.1 

20.549 

17.648 

2.938 

1.704 


-1.336 

710 

05 

31 

23.50 

+34 

14 

02.0 

21.590 

18.543 

3.048 



-1.852 

711 

05 

31 

23.51 

+34 

18 

16.3 

20.080 

18.087 

1.985 

1.519 


0.055 

712 

05 

31 

23.53 

+34 

14 

44.8 

21.151 

18.865 

2.286 




713 

05 

31 

23.53 

+34 

13 

41.5 

17.125 

15.433 

1.629 

1.346 

0.648 

0.023 

714 

05 

31 

23.55 

+34 

15 

05.5 

16.541 

15.444 

1.091 

0.782 

0.479 

0.313 

715 

05 

31 

23.57 

+34 

13 

09.2 

21.865 

18.669 

3.196 



-1.347 

716 

05 

31 

23.57 

+34 

13 

35.3 

19.420 

17.059 

2.377 

1.811 



717 

05 

31 

23.60 

+34 

10 

59.6 

20.964 

17.982 

2.958 

1.642 


-0.206 

718 

05 

31 

23.60 

+34 

13 

54.2 

19.331 

17.506 

1.826 

1.543 


-0.105 

719 

05 

31 

23.60 

+34 

14 

10.2 

21.592 

18.259 

3.333 




720 

05 

31 

23.61 

+34 

10 

32.2 

21.686 

18.904 

2.787 




721 

05 

31 

23.64 

+34 

11 

12.7 

21.413 

19.373 

2.021 




722 

05 

31 

23.68 

+34 

10 

40.7 

15.000 

14.352 

0.640 

0.472 

0.201 

0.314 


ha for NGC 1931 


€V 


Bn 

^B-V 

^U-B 

^H-C 


^obs 



2MASSID 

Ha 

0.007 

0.031 

0.058 



0.079 

2 

3 

2 

0 

0 

1 

05312347+3412276 

H 

0.041 

0.037 

0.071 

0.221 


0.099 

2 

3 

2 

2 

0 

1 

05312345+3410471 

H 

0.001 

0.013 

0.013 

0.062 

0.083 

0.073 

3 

4 

3 

2 

2 

1 

05312347+3413281 

- 

0.085 

0.029 

0.050 




2 

2 

2 

0 

0 

0 

05312349+3412182 

- 

0.001 

0.045 

0.056 

0.105 


0.058 

2 

4 

2 

2 

0 

1 

05312347+3413060 

H 

0.084 

0.054 

0.100 



0.109 

1 

1 

1 

0 

0 

1 

- 

- 

0.027 

0.001 

0.021 

0.062 


0.058 

2 

3 

2 

2 

0 

1 

05312351+3418161 

- 

0.104 

0.062 

0.121 




1 

1 

1 

0 

0 

0 

05312356+3414448 

- 

0.020 

0.019 

0.004 

0.031 

0.014 

0.013 

4 

4 

4 

4 

4 

2 

05312353+3413415 

- 

0.008 

0.002 

0.010 

0.006 

0.009 

0.012 

4 

4 

4 

4 

4 

2 

05312355+3415056 

- 

0.116 

0.047 

0.125 



0.154 

1 

1 

1 

0 

0 

1 

05312357+3413097 

- 

0.045 

0.004 

0.044 

0.057 



4 

4 

4 

2 

0 

0 

05312356+3413353 

- 

0.047 

0.019 

0.023 

0.131 


0.102 

2 

4 

2 

1 

0 

1 

- 

- 

0.001 

0.038 

0.044 

0.100 


0.072 

4 

4 

4 

2 

0 

1 


- 

0.104 

0.031 

0.109 




1 

2 

1 

0 

0 

0 

O 

- 

0.017 

0.013 

0.031 




2 

2 

2 

0 

0 

0 

05312361+3410320 1 

- 

0.041 

0.068 

0.119 




2 

2 

2 

0 

0 

0 

- 

- 

0.005 

0.001 

0.004 

0.003 

0.005 

0.007 

4 

4 

4 

4 

4 

2 

05312367+3410407 

- 
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A. TRANSFORMATION RELATIONS FOR THE MONT4K CCD 
CAMERA OF THE KUIPER 61” TELESCOPE 


The Mont4k photometric system of the Kuiper 61” telescope at SO was used in the 
SOS project for the hrst time. It is necessary to obtain reliable transformation relations of 
the instrumental system to produce homogeneous standardized photometric data for many 
open clusters. In this appendix we describe the derivation of the transformation relations 
using a large number of the standard stars observed in 2011 - 2013 (281 in Bessell U, 119 
in SDSS u', 402 in B, 436 in V, and 530 in I). The SDSS u' hlter was used on several 
nights in 2012, instead of the Bessell U hlter, and so the transformation relation was also 
addressed for the forthcoming papers. 


The Mont4k tasks installed in the observing system automatically corrected for 
variations in the bias level and crosstalk for individual images, converting from the raw 
extended FITS into normal FITS images. In order to remove remaining instrumental 
artifacts we obtained a number of bias, sky hat, and shutter shading images during 
observing runs. The bias level was very small (a few ADU) because the overscan correction 
had already been done. After subtracting the tiny residual bias, all target images were 
divided by the master hat images in the same passbands. The Mont4k CCD camera 
system is equipped with an iris-type shutter. This type of shutter imprints an iris pattern 
on short-exposure images. In addition, there could be either a systematic gain or loss in 
the observed hux depending on the speed of the shutter. These behaviours are crucial 
for obtaining reliable magnitudes in short exposure images. We carefully investigated the 
shutter shading with various exposure tin ie f0.5s. Is , 2s, 3 s, 4s, 5s, 7s, 10s, 15s, and 20s) 


according to the procedures delineated by 


Lim et ah 


(120081) and found that images exposed 


for less than 7s showed the iris pattern. Images with exposure times less than 2s revealed 
a systematic gain variation of more than 0.5 percent (up to 4.6 percent at 0.5s) in the 
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observed flux. It implies that the shutter operates more slowly for exposure time shorter 
than 2s. Based on these results we corrected the shutter shading effect on short exposure 
images (0.5 - 7s). 

We performed aperture photometry for the standard stars with an aperture size of 14.0 
arcsec (16.3 pixels). The enormous amount of photometric data obtained over 11 nights 
in 2011 - 2013 allowed us to obtain reliable transformation relations. The instrumental 
magnitudes can be transformed to the standard magnitudes and colors using the following 
equation: 


Mx = mx-{kix-k2xCo)-X+7]x-Co+ax-UT+(3x-XccD+l\-^ccD+^\-yccD+(^-yccD+C,\ (Al) 


where /^a, 7a, 5a, and ex represent the coefficients of spatial variation in magnitude. Other 
symbols are the same as in Equation 1. The Xccd and Yccd coordinates were normalized 
by 1,000 pixels [xccd = 000, yccD = hccD/1,000). Figure lA.II shows the spatial 

variations in the measured magnitudes with respect to Xccd and Yccd coordinates for 
each passband. Since there are five uncountable columns in the center of the normal FITS 
images, the Xccd coordinate of stars away from the center (682*^ pixel) in a positive 
direction was increased by 5 pixels. The form of the spatial variations can be represented 
by a second-order polynomial. T he pattern may relate 


the 


20i; 


<uiper 61” telescope (see also 


Magnier fc Cnillandre 


to the 

shape of the 

foca 

plane of 

to 

o 

o 

_ 

Sung et al. 

2008 

him et al. 


). In most cases, the amount of the variation from the CCD center was less than 0.03 


mag. The maximum difference was up to 0.05 mag in the Bessell U band. We performed 
the second-order polynomial fitting to the data, presenting the coefficients of the spatial 
variation in Table lA.ll In addition, a linear variation along the Yccd coordinate in the I 
band was found between the 0 - 500*^ pixels, but there was no remarkable spatial variation 
on the other side. 
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The atmospheric extinction and transformation coefficients were determined using a 
weighted least-square method after correcting for the spatial variations. Figure |X]2] shows 
the transformation relations for the Mont4k CCD camera system. A combination of a few 
straight lines constitutes the transformation relations. Since there was only one extremely 
red star (1/ — J > 3 ) in the I transformation with respect toV — I color, the transformation 
relation may be uncertain in that color range. The V transformation with respect to V — I 
shows a plateau for the reddest stars {V — I > 1.55). The TiO band spectral features 


in late-tvpe stars 

are related to t 

lis pi 

2000; 

him et ah 

2009; 

him et ah 

2013) 


ateau as seen in other studies flSung fc Bessell 


20131) . The B transformation was approximated by a 


combination of two straight lines, and the coefficients are very small. The transformation 
relations of the Bessell U and SDSS u' hlters exhibit opposite trends t o each other. A 


non-linear 


jehaviour was fou nd in the U transformation. According to 


fl2nnnh and 


him et ah 


Sung fc Bessell 


(120091) this aspect is caused by the Balmer discontinuity in the fluxes 


of B-F-type stars and difference between the response functions of the Johnson-Cousins 
standard system and natural photometric systems. Such a non-linearity can be corrected 
using dereddened B — V colors as shown in the lower-right panels of Figure lA^ The aspect 
of the non-linear correction term appears as an opposite trend in the Bessell U and SDSS 
u' transformation relations. The transformation coefficients with respect to given colors are 
summarized in Table lAAl 

On the other hand, we found that standard stars with an intrinsically red color 
appeared systematically brighter than other stars in the Bessell U hlter after correcting 
with all the coefficients above (see the upper panel in Figure lAT^ . Unlike the Bessell 
U hlter used at MAO, the U hlter attached to the Mont4k CCD camera exhibits the 
so-called red leak. Since the contribution of red light drastically increases for U — B > 1, 
a careful correction for the red leak ehect was carried out with respect to the V — I color 
as shown in the upper panel of Figure IA731 The transformation of the SDSS u' hlter also 
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required an additional correction term for extremely red stars. Because the transmission 
function of the SDSS u' filter was shifted toward shorter wavelength to avoid the Balmer 
discontinuity, one can expect that the u' magnitude of late-type stars may be fainter than 
the Johnson U magnitude. Thus, a negative correction term for red stars {V — I > 1.5) was 
required as shown in the lower panel of Figure EHl We note that all the coefficients in the 
transformation equation were obtained from iteration of the procedures described above. 
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0.05 
AI 0 

-0.05 
0.05 
AV 0 

-0.05 
0.05 
AB 0 

-0.05 
0.05 
AU 0 

-0.05 
0.05 
Au 0 

-0.05 

0 10 1 
Xccd/1000 Y,V1000 

Fig. A.l.— Spatial variation of the measured magnitudes with respect to the Xccd (left) 
or Yccd (right) coordinates of the Mont4k CCD chip. The Xccd coordinates of stars in a 
positive direction from the CCD center {Xccd = 682) was increased by 5 pixels because of 
five uncountable columns in the center of the normal FITS images. All the pixel coordinates 
of the standard stars on the CCD chip have been divided by 1,000. Most spatial variations 
appear in the form of a second-order polynomial equation. The variation of the I magnitude 
along the Yccd can be approximated by a combination of two straight lines. The maximum 
difference between the center and edge of the CCD is less than 0.05 mag. 
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Fig. A.2.— Transformation relations of the Mont4k CCD camera system with respect to rel¬ 
evant colors. The size of circles is proportional to the nnmber of observations and the degree 
of precision of the standard stars. The error bar is the standard deviation of the observed 
magnitnde. The combinations of straight lines in each panel represent the transformation 
relation. The lower-right panels show a non-linear correction term against dereddened B — V 
color in the U band. See the main text for details. 
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0 1.5 3 

V-I 


Fig. A.3.— Additional correction terms in the U band for the intrinsic red stars {U — B > 
1.0). Upper: Red leak effect found in the Bessell U hlter attached to the Mont4k CCD 
camera system. The contribution of the red light can be removed by a linear regression for 
stars (y — I > 1.79). Lower: More flux in the U band is required for late-type stars because 
our SDSS u' hlter was designed to avoid the Balmer discontinuity toward shorter wavelength. 
An additional correction can be made with a straight line for red stars {V — I > 1.5). 
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Table A.l. 


Spatial Variation Coefficients of the Mont4k CCD Camera 


Filter 

/3 

7 

& 

e 

Remark 

I 

0.062 

-0.086 


-0.052 

Yccd < 500 pixels 

V 

0.024 

-0.029 

0.046 

-0.066 


B 

0.029 

-0.029 

0.043 

-0.063 


U 

0.042 

-0.052 

0.058 

-0.097 


u' 

0.030 

-0.049 

0.047 

-0.066 
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Table A.2. Transformation Coefficients of the Mont4k CCD Camera 


Filter 



Color range 




V\,l 

Vx ,2 

V\,3 

V\,4 

Correction term 

I 

V-I <3 

V-I>3 





0.015 ±0.001 

0.066 




V 

V - I <0.2 

0.2 < F - / < 0.41 

0.41 <V -I < 1.55 

V -I > 1.55 



-0.080 ±0.010 

0.000 

-0.054 ±0.002 

0.000 



B -V <0.3 

0.3 < B - V < 0.45 

V - I > 0.45 




-0.077 ± 0.006 

0.000 

-0.047 ±0.001 




V - R< 0.13 

0.13 <V - R< 0.25 

0.25 <V - R< 0.45 

V -R>0.9 



-0.150 ±0.013 

0.000 

-0.094 ± 0.003 

0.000 


B 

B-V <0.5 

B -V >0.5 





0.013 ±0.002 

-0.005 ±0.001 




U 

U - B < -0.57 

-0.57 <U - B < 0.67 

U - B> 0.67 


V -I > 1.79 


0.253 ±0.016 

0.000 

0.069 ± 0.007 


0.123 ± 0.019 

u' 

U -B <0. 

0. <U - B <0.A 

U - B >0.4 


V -I>1.5 


-0.179 ± 0.004 

0.043 

-0.111 ±0.004 


-0.120 ±0.007 










